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C3 ABSTRACT 

The extragalactic background light (EBL) is of fundamental importance both for 

_a_ understanding the entire process of galaxy evolution and for 7-ray astronomy, but 

£> the overall spectrum of the EBL between 0.1 - 1000 /zm has never been determined 

^^ directly from galaxy spectral energy distribution (SED) observations over a wide 

l/~j redshift range. The evolving, overall spectrum of the EBL is derived here utilizing 

a novel method based on observations only. This is achieved from the observed 
t—I evolution of the rest -frame A'-band galaxy luminosity function up to redshift 4 

1^ ( |Cirasuolo et al.|2010 ), combined with a determination of galaxy SED-type fractions. 

f^ These are based on fitting SWIRE templates to a multiwavelength sample of about 

f*^ 6000 galaxies in the redshift range from 0.2 to 1 from the All-wavelength Extended 

y— I Groth Strip International Survey (AEGIS). The changing fractions of quiescent 

galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) 
galaxies in that redshift range are estimated, and two alternative extrapolations of 
SED-types to higher redshifts are considered. This allows calculation of the evolution 
^ of the luminosity densities from the UV to the IR, the evolving star formation rate 

density of the universe, the evolving contribution to the bolometric EBL from the 
different galaxy populations including AGN galaxies and the buildup of the EBL. 
Our EBL calculations are compared with those from a semi-analytic model, from 
another observationally-based model and observational data. The EBL uncertainties 
in our modeling based directly on the data are quantified, and their consequences 
for attenuation of very high energy 7-rays due to pair production on the EBL are 
discussed. It is concluded that the EBL is well constrained from the UV to the 
mid-IR, but independent efforts from infrared and 7-ray astronomy are needed in 
order to reduce the uncertainties in the far-IR. 

Key words: galaxies: evolution - galaxies: formation - cosmology: observations - 
diffuse radiation - infrared: diffuse background 
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1 INTRODUCTION 

The formation and evolution of galaxies in the universe is ac- 
companied unavoidably by the emission of radiation. All this 
radiated energy is still streaming through the universe, al- 
though much is now at longer wavelengths due to redshifting 
and absorption/re-emission by dust. The photons mostly lie 
in the range ~ 0.1-1000 /mi, i.e., ultraviolet (UV), optical, 
and infrared (IR), and produce the second most energetic 
diffuse background after the Cosmic Microwave Background, 
thus being essential for understanding the full energy bal- 
ance of the universe. We will account in this work for the 
radiation accumulated by star formation processes through 
most of the life of the universe, plus a contribution from ac- 
tive galactic nuclei (AGNs) to this wavelength range, known 
as the diffuse extragalactic background light (EBL). 

The direct measurement of the EBL is a very difficult 
task subject to high uncertainties. This is mainly due to 
the contribution of zodiacal light, some orders of magnitude 
larger than th e EBL (e.g., |Hauser fc Dwek||200l| |Chary| 
& Pope||2010|. There are some measurements in the opti- 
cal ([B ernstein 2007[) and in the near-IR (e.g., Cambresy et 
|al.||200 1; Mats umoto et aL]|2005[ ), but there is not general 
agreement about the reliability of these data sets (Mattila 



20061. In addition, these near-IR data appear to give inten- 



sity levels for the EBL in contradiction with the observa- 
tion of very high energy ( VHE, 30 GeV-30 TeV ) photons 
from extragalactic sources ( Aharonian et al.|2006 Mazin & 
|Raue||2007] | Albert et aT1|2008[ ). Little is known about the 
mid-IR from direct detection due to the higher contamina- 
tion from zodiacal light at those wavelengths. Measurements 
with the Far-Infrared Absolute Spectrometer (FIRAS) in- 
strument on board the Cosmic Background Explorer, in the 
far-IR ( Hauser et al.|1998 Lagache et al.|2000 l, are thought 
to be more reliable. Other observational approaches set re- 
liable lower limits on the EBL, such as measuring the inte- 



grated light from discrete extragalactic sources (e.g., Madau 
|fc Pozzetti||2000| [Fazio et al.|2004"| ). 

There are also other authors that focus on studying 
galaxy properties based on EBL results ( Fardal et al.|2 007), 
or on modeling a region of the EBL spectrum ( Younger & 



Hopkins 2010). On the other hand, there are phenomeno- 



logical approaches in the literature that predict an overall 
EBL model (i.e., between 0.1-1000 /mi and for any redshift). 
These are basically of four kinds: 

(i) Forward evolution, which begins with cosmological ini- 
tial conditions and follows a forward evolution with time by 
means of semi-analytical models (SAMs) of galaxy forma- 
tion, e.g., Primack et al.| (19991, Somerville et al. (20101 



(hereafter, SGPD10) and |Gilmore et al.| < |2010| > (hereafter, 
GSPD10). 

(ii) Backward evolution, which begins with existing 
galaxy populations and extrapolates them backwards in 
time, e.g., 



Malkan & Stecker (19981, Stecker, Malkan & 



Scully (20061, Franceschini, Rodighiero fc Vaccari ( 2008 1 



(hereafter, FRV08). 

(iii) Evolution of the galaxy populations that is inferred 
over a range of redshifts. The galaxy evolution is inferred 
here using some quantity derived from observations such 
as the star formation rate (SFR) density of the universe, 
e.g., Kneiske et al. (2002), Finke, Razzaque & Dernier 



(iv) Evolution of the galaxy populations that is directly 
observed over the range of redshifts that contribute signifi- 
cantly to the EBL. The present paper, which we term em- 
pirical, belongs in this category. 

The type (i) SGPD10 and GSPD10 models discuss 
the same galaxy formation SAM but in different contexts: 
SGPD10 contains details of the model used in calculating the 
bolometric luminosity history of the universe and compari- 
son with data, and GSPD10 focuses on the derived EBL and 
7-ray attenuation. The SGPD10-GSPD10 model is based 
on an updated version of the semi-analytic theoretical ap- 
proach described in Somerville et al. ( 2008 1 from the growth 



of super-massive black holes and their host galaxies within 
the context of the hierarchical Lambda Cold Dark Matter 
(ACDM) cosmological framework. This is based in part on 



Somerville & Primack (19991, Somerville, Primack, & Faber 
(|2001|), and in the simulations summarized by Hopkins et al. 



2008a), and Hopkins et al. (2008b). We consider that these 



(20101, Kneiske & Dole (20101 



types of models are complementary to the observational ap- 
proach taken here. 

We consider the type (ii) FRV08 model the most com- 
plete observationally-based work of those mentioned above. 
They base their EBL modeling on galaxy luminosity func- 
tions (LFs), quantities which are directly observed and well 
understood. FRV08 exploit a variety of data to build evolu- 
tionary schemes according to galaxy morphology. They ac- 
count for the contribution from early, late-type galaxies and 
a starburst population to the EBL. They use observed near- 
IR LFs from the local universe to z = 1.4 for describing the 
early and late-type galaxies. For the starburst population 
they use an optical and only local LF. Different prescrip- 
tions are used to extrapolate the evolution of the different 
morphological types to higher redshifts, and corrections to 
fit their results to other observational data are applied. 

Type (iii) models are not directly based on galaxy data. 
Instead, they are built from some parametrization of the 
history of the SFR density. This is a quantity derived us- 
ing several different methods, each of which have different 
and signitficant uncertainties and biases. The SFR density 
is combined with uncertain assumptions about the emitted 
galaxy spectral energy distribution (SED) evolution as well. 

Our type (iv) EBL estimates (the first approach in this 
category) will be compared in detail with the type (i) for- 
ward evolution semi-analytical galaxy formation model by 
SGPD10 and GSPD10, and the type (ii) oservationally mo- 
tivated model by FRV08. The other works mentioned are 
briefly compared with our EBL calculations in Sec. [6] 

Our aim in this paper is to develop an EBL model that 
is as observationally-based and realistic as possible, yet fully 
reproducible, including a quantitative study of the main un- 
certainties in the modeling that are directly due to the data. 
This constrains the range of the background intensity and its 
implications to 7-ray astronomy. One important application 
of the EBL for 7-ray astronomy is to recover the unattenu- 
ated spectra of extragalactic sources. Our goal is to measure 
the EBL with enough precision that the uncertainties due to 
the EBL modeling, in these recovered unattenuated spectra, 
are small compared with other effects such as systematic un- 
certainties in the 7-ray observations. Examples of this are 
discussed in Sec. [5] 

Our model is based on the rest-frame A'-band galaxy LF 
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Figure 1. Histogram of the number of galaxies versus redshift of 
our sample in the four redshift bins considered in this work. 
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Table 1. The photometric bands in our galaxy sample. For each 
we show the effective wavelength, the data source, the require- 
ment for that band to be included for a given galaxy in our sam- 
ple (det: a detection in this band is required; obs: observation in 
this band is required, but not necessarily a detection; ext: this 
band is considered extra information when available), and the 5<r 
upper limit in that band in cases where there is no detection. 



by Cirasuolo et al. ( 2010 1 (hereafter, C10) and on multiwave- 



length galaxy data from the All-wavelength Exte nded Groth 
Strip Internatio nal Survey (AEGI971 |Davis et al.|2007||New- 
man et al.|2010 1 of about 6000 galaxies in the redshift range 
0.2-1. These data sets are put together in a very transpar- 
ent and consistent framework. The C10 LF is used to count 
galaxies (and therefore to normalize the total EBL inten- 
sity) at each redshift. The LF as well as our galaxy sample 
are divided in three magnitude bins according to the abso- 
lute rest-frame A-band magnitude i.e., faint, middle, and 
bright (defined quantitatively later). Within every magni- 
tude bin a SED-type is statistically attached to each galaxy 
in the LF assuming SED-type fractions that are function of 
redshift within those magnitude bins. This is estimated from 
fitting our AEGIS galaxy sample to the 25 galaxy SED tem- 
plates from the SWIRE library. Then, luminosity densities 
are calculated from these magnitude bins from every galaxy 
population at all wavelengths, and finally all the light at all 
redshifts is added up to get the overall EBL spectrum. The 
results are linked with 7-ray astronomy and with the current 
understanding on galaxy evolution. 

The paper is organized as follows: Section [2] describes 
the LF, our multiwavelength galaxy catalogue and the 
galaxy templates. Section 13] explains our methodology. The 
results for galaxy SED-type fractions, luminosity densities, 
SFR densities, EBL buildup, and EBL intensities are given 
in SectionH] Section[5]shows the attenuation computed from 
our EBL model for some VHE sources taken from the liter- 
ature. In Section [6] our results are discussed including a de- 
tailed study on the uncertainties from the modeling, a com- 
parison between our observationally-based EBL and that 
given by theoretical SAMs of galaxy formation. Finally, in 
Section [7] a summary with our main results and conclusions 
is presented. 

Throughout this paper, a standard ACDM cosmol- 
ogy is assumed, with matter density Q m — 0.3, vacuum 
energy density S1a = 0.7, and Hubble constant Hq — 
70 km s _1 Mpc _1 . 



1 http://aegis.ucolick.org/ 

© 2010 RAS, MNRAS 000,[l}J25] 



2 DATA DESCRIPTION 

2.1 A'-band galaxy luminosity function 

The evolving galaxy LF in rest-frame A'-band provided by 
C10 from z = to 4 is used. This evolving LF is the most 
accurate measurement to date of cosmological galaxy evo- 
lution in the near-IR, where dust absorption is less severe 
than in optical bands. The fc-corrections in this band are 
less severe than in the optical as well. The choice of the C10 
LF to normalize the model is also based on the smooth and 
well-studied shape of the galaxy SEDs in the near-IR, unlike 
others in UV or mid-IR wavelengths. 

The resulting evolving LF is based on the UKIDSS Ul- 
tra Deep Survey ( |Lawrence et al.||2007[ ), which has a large 
area and depth, and hence reduces the uncertainties due 
to cosmic variance and survey incompleteness. We refer the 
reader interested in details to that work. It is important 
to note that they give a parametrization of the evolution 
of the LF corrected from incompleteness and fitted by a 
Schechter function ( |Schechter|1976| over redshift, <I>(Mjj-, z), 
where Mjj- is rest-frame A-band absolute magnitude at red- 
shift z. The strongest assumption that they make is to keep 
constant the faint-end slope a in their parametrization. 



2.2 Galaxy sample description 

A multiwavelength galaxy catalogue built from AEGIS 
( Davis et al.|2007 Newman et al.|2010 1 for this work is used. 
This catalogue contains 5986 galaxies, all in the Extended 
Groth Strip (EGS). It is required that every galaxy in the 
sample have 5<r detections in the B, R, I, Ks and IRAC 1 
bands, and observations (but not necessarily detections) in 
the IRAC 2, 3, 4 and MIPS 24 bands (see Table [TJ. These 
5a upper limits are given by the following fluxes: 1.2, 6.3, 
6.9 and 30 fj, Jy, for IRAC 2 3, 4 a nd MIPS 24 respectively, 
according to|Barmby et al.| ((2008) for the IRAC bands and 
Dickinson et al.| (|2007[ ) for MIPS 24. They are also summa- 
rized in Table 111 In addition, 1129 of these galaxies have 
GALEX detections in the far-UV and 2345 galaxies in the 
near-UV. In our sample, 4376 galaxies have the highest qual- 
ity spectroscopic redshifts measured by the Deep Evolution- 
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Figure 3. Colour-magnitude diagram in the same four different redshift bins showing the galaxies of our sample after the cuts explained 
in Scc. |4.l| for the magnitude bins defined in the text for the integrals in Eq. [3] It is over plotted the luminosity function by |Cirasuolo| 
|et al.| ( |2010^ in the mean of every redshift bin with arbitrary units in the logarithmic j/-axis. The colour code is the same that in Fig. |2| 
Magnitudes are in AB magnitude system. 
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ary Exploratory Probe 2 team (DEEP2 DR3, Newman et al 



2010]), with the DEIMOS spectrograph ( |Faber et al.||2003 1 
on the Keck II telescope in an area of about 0.7 deg in 
the sky. All the other galaxies in the sample (1610 galax- 
ies) have secure photometric redshifts, more than 80% with 
uncertainty in redshift less than 0.1. The redshift covered is 
between 0.2-1 (almost 60% of the age of the universe) for a 
total sample of 5986 galaxies. For our purpose we will not 
distinguish between spectroscopic or photometric redshifts. 
This assumption will be discussed in Sec. |6.1| 

The optical photometry (B, R and I bands) was taken 
from imaging with the CFH12K camera ( |Cuillandre et| 
al. 2001 1 on the Canada-France-Hawaii Telescope (CFHT) 
3.6 m. The integration time for these observations was 1 hr 
in B and R, and 2 hr in I, per pointing. More details may 



be found in Coil et al. ( 2004 1 



The near-IR photometry in the Ks band is from the 
Wide-field Infrared Camera (WIRC, | Wilson et aL|[2003] | 
camera on the Hale 5 m telescope at the Palomar obser- 
vatory. This data set is the most restrictive constraint on 
the area of our sample, therefore our galaxy catalogue is Ks 
limited. The ECS field was surveyed to different depth for 
different sub-regions up to Ks — 22.5 in the AB magnitude 
system. The details may be found in |Conselice et aLJ ([2008). 

The mid-IR data are from the IRAC and MIPS cameras 
on board the Spitzer Space Telescope. The details may be 
found in |Barmby et al.| ( |2"008| ) and in |Dickinson et al.|f2007| ) 
describing the FIDEL survey, the source of our 24 /xm data. 

In addition, some data in the UV in two different bands 
0.1530, 0.2310 /im from the Galaxy Evolution Explorer 



(GALEX, Morrissey et al. 20071 are included in our cat- 



alogue. This data set is part of the GALEX Deep Imaging 
Survey and the details may be found in |Salim et al.| f2009). 

Source catalogue from each of these imaging data sets 
where cross-matched using a Bayesian method, which took 
into account prior information from the surface densities of 
sources in each band ( |Huang et al.|2010 |. The IRAC 1 data 
were used as the primary reference catalogue. 

It is important to note that all our data are public, 
except the MIPS 24 photometry, the cross-match catalogue, 
and the photometric-redshift catalogue ( Huang et al.| 2010). 
These will be released to the public soon. 

The histogram of the redshift distribution of the AEGIS 
sample is shown in Fig.fllin the four redshift bins considered 
in our calculations. Note the larger number around z ~ 0.7, 
mainly due to the weighting scheme of the DEEP2 survey, 
which tends to select galaxies at z > 0.7 based on colour- 
colour criteria, plus the effect in the opposite direction of 
losing faint galaxies at higher redshifts. 

In order to calculate the absolute magnitudes in U, B 
and if-band, we have computed the best-fitting template 



taken from the Bruzual & Chariot ( 2003 1 stellar population 
models to the data in the photometric bands B, R, I, Ks, 
IRAC 1 and FUV, NUV and IRAC 2 when available, using 
the code FAST (Fitting and Assessment of Synthetic Tem- 



plates, see the Appendix in Kriek et al. 2009 f or details). 
FAST makes a \ 2 minimization from a grid of Bruzual & 



Chariot (20031 models. We chose a stellar initial mass func- 



tion (IMF) given by Chabrier ( 2003J , an exponentially de- 
clining SFR ~ exp(— t/r) with r ranging from 10 7 -10 10 Gyr, 
(the same range for the ages) , metallicities by mass fraction 
in the range 0.004-0.050 (solar metallicity is 0.02 in these 
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Figure 4. Contribution from the three different magnitude bins 
defined in Sec.|3]to the total of the co-moving rest-frame K-hand 
luminosity density calculated directly from the luminosity func- 
tion (LF) by |Cirasuolo et alT] ( |2010[ | . The bulk of the light comes 
from the middle and bright-end of the LF, where the Schechter 
parameter L* is. Note the increment with redshift of the bright- 
end contribution which decreases the impact of a possible color- 
selection effect or mis- typing (see Sec. |6.1y at the highest redshift 
in our galaxy sample. 



units), and optical extinction A v from 0-10 following the 
|Calzetti et al.| ( |2000[ ) extinction law. We calculate the ab- 
solute magnitudes from the best-fitting model using the U 
Bessel filter, the B filter from CFHT12K, and the same K- 
band filter from the UKIDSS survey, the same filter where 
the LF from C10 was estimated. All the transmission curves 
for these filters can be found in the default distribution of 
Le PHARE. 

The sample was not corrected for incompleteness. How- 
ever it is estimated here how this affects our results. The 
colour-dependent incompleteness of the DEEP2 survey was 



studied in Willmer et al. (20061. They estimated a relation 



between the rest-frame U — B colour versus the absolute 
magnitude in the B-band Mb for which galaxies fainter 
than this relation have colour-dependent incompleteness. We 
show in Fig. [2] colour-magnitude diagrams of our AEGIS 
galaxy sample for four different redshift bins. The black line 
is taken from Fig. 4 in Willmer et al. (20061. Galaxies lo- 



cated to the right of this line are likely missing. This figure is 
colour coded according to the calculated best-fitted SWIRE 



template (see Sec. 2.3 and Sec. 4.1l. The number of galax- 



ies lying to the right of the relation, thus suffering colour- 
dependent incompleteness, are only of 1.8%, 2.3%, 7.3% and 
9.3% for the different redshift bins presented in Fig. [2] 

Fig. [5] shows rest-frame U — B versus absolute magni- 
tude in the /("-band in the three magnitude bins considered 
in this work to show an estimation on the galaxy number in 
each bin and their SED-types. We will describe this figure 
in the context of cosmological evolution in Sec. |6.4| 

Thus to recap, the normalizations of the EBL in our 
model is given by the i("-band LF of C10, and our galaxy 
SED-type fractions give the relative contribution of every 
galaxy type to the total luminosity density. The assignment 
of SED-types to the galaxy population at a given redshift is 
done individually for three ranges in K-band absolute mag- 
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nitude, as it will be discussed in Sec. [3] Moreover, most of 
the contribution to the EBL (between 70-90%) conies from 
around the knee of the LF (__■* according to the Schechter 
parametrization) as shown in Fig. El for the rest-frame K- 
band luminosity density (calculated directly from the inte- 
gration of the CIO LF) and not from the faintest galaxy 
population where we suffer some small colour-dependent in- 
completeness. Fig HI also shows that the contribution from 
the bright-end increases with redshift decreasing the impact 
of any colour-dependent effect. As the remaining colour- 
independent incompleteness does not have any effect on the 
galaxy SED-type fractions in our model or the overall nor- 
malization (which is set by the if-band LF), we conclude 
that our results are quite robust and the effect from incom- 
pleteness in our sample is minimal. 



2.3 Galaxy spectral energy distribution library 

The galaxy SEDs fo und in the SWIRE template librar^ 
( Polletta et al.|20 07) are used. This library contains 25 tem- 
plates, representative of the local galaxy population, defined 
as 3 elliptical galaxies, 7 spirals galaxies, 6 starbursts, 7 
AGN galaxies (3 type I AGN, 4 type II AGN), and 2 com- 
posite (starburst+AGN) templates all covering the ~ 0.1- 
1000 )im wavelength range. The elliptical (quiescent), spi- 
ral (star-forming) and starburst (very star-forming) IR tem- 
plates were generated with the GRASIL code ( |Silva et al.| 
1998) based on observations. The 7 spirals range from early 
SO- 



to late types (i.e., SO - Sdm). The starburst templates corre- 
spond to the SEDs of NGC 6090, NGC 6240, M 82, Arp 220, 
IRAS 22491-1808, and IRAS 20551-4250. In all of the spi- 
ral and starburst templates, the spectral region between 5- 
12 /j,m, where many broad emission and absorption features 
are observed, was replaced by observed IR spectra from the 
PHT-S spectrometer on board the Infrared Space Obser- 
vatory and from IRS on Spitzer. Some examples of these 
templates are shown in Fig. [5] 

We are aware that these templates do not include high- 
redshift galaxies (z > 0.3). This effect will be taken into 
account in a future version of our EBL model when higher- 
redshift galaxy SEDs are released, including very vigorous 
starburts and AGNs not present in the local universe. The 
limitation of using local templates for high-redshift galax- 
ies has been addressed in some works such as Mu rphy et| 
al. (20091, where they conclude that for ultra-luminous IR 



galaxies (ULIRGs) in the redshift range between 1.4 < z < 
2.6, IR luminosities are overpredicted when they are derived 
only using MIPS 24 photometry, thus showing a different 
behaviour than local ULIRGs. It has been also shown by 



Menendez-Delmestre et al. ( 2009 1 that submillimeter galax- 
ies in the redshift range between 0.65 < z < 3.2 show 
a larger polycyclic-aromatic-hydrocarbon (PAH) emission 
than local analogs, suggesting a more extended dusty star- 
forming region than seen in local ULIRGS. 
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Figure 5. Spectral energy distributions for some galaxy tem- 
plates from the SWIRE library. We show here (from the bottom 
to the top) an early type quiescent galaxy (EU13), a very late 
star-forming galaxy (Spi4), a starburst galaxy (122491) and two 
different AGN galaxies: a Seyfert II, and a quasi-stellar object 
type I (QSOl). The j/-axis is in arbitrary units. 



3 METHODOLOGY 

The empirical approach of the EBL evolution directly ob- 
served over the range of redshifts that contribute signifi- 
cantly to the EBL is followed. This is type (iv) according 
to the classification given in Sec. [I] As briefly explained in 
Sec. [I] our aim is to calculate the EBL integrating over red- 
shift luminosity densities. These quantities are estimated at- 
taching statistically SEDs to the galaxies given by the LF 
by C10 in three different magnitude bins. This is achieved 
using galaxy SED-type fractions between z = 0.2 — 1 by 
finding the best-fitting template of the 25 SED templates in 
the SWIRE library describing every galaxy in the AEGIS 
galaxy sample. Two different extrapolations for the galaxy 
SED-type fractions for z > 1 are assumed leading to the 
same evolving EBL intensity from the UV to the mid-IR 
but different far-IR. 

The Le PHARE v2.2 (Photometric Analysis for Red- 
shift Estimations) code is used to find the best-fitting 
SWIRE SED template for every gala xy in the sample. Le 

( |Arnouts et al.||2002 



PHARE is a publicly available coda 3 

|Ilbert et al.||2006"| ) mainly aimed to calculate photometric 
redshifts, but with the possibility to find the best-fitting 
template (among any library introduced as input) for galax- 
ies with known redshift. Le PHARE makes use of a \ 2 fitting 
procedure weighted from normalizations in every detected 
bands, and with the possibility to set upper limits for fluxes 
in some bands based on non-detections. From the fact that 
we have required observations in several bands to build our 
catalogue, we set for every galaxy in the fitting procedure 
upper limits on the bands where there is no a 5<r detection. 
The information at all bands is used in the fitting. 

For every galaxy, templates are rejected if they predict 
a flux in a given band that is higher than the upper limit 



2 http://www.iasf- milano.inaf.it/'~polletta/templates/ 
swire_templates.html 



3 http://www.cfht.hawaii.edu/~arnouts/LEPHARE/cfht_lephare/ 
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1 io° K^icr 1 io° ic^ict 1 io° ioHct 1 io° io 1 
Xrest [/inn] \ rest [/iin] \ rest [/iin] \ rest [/iin] 



Figure 6. Examples of the best fits (upper panel), fits around X 2 . e d = ^® (second upper panel), fits around X^ ed =30 (second lower panel) 
and the worst fits (lower panel). The columns are from left to right: quiescent, star-forming galaxies, starbursts, and AGN galaxies. The 
AEGIS identification number is shown for the galaxy along with X^-ed & ven by the fitting code Le PHARE described in Sec. 3 The 
information at all bands is used in the fitting. 



for that band. The equations used for the fitting procedure 
are shown in Eq. IT] with the parameter s given by Eq. [2] 



2 
X = 



E 



Fob 



SJ^temp,i 



E 



■F obs.j -Fte 



•-)' 






F, 



(1) 



(2) 



The notation is the following: i refers to a given band, 
j to the band used for the scaling, ¥ a bs is the observed flux, 
Ftemp is the flux from the templates, a is the la statisti- 
cal uncertainty of the photometric measurement, and s is 
the scaling factor that is chosen to minimize the x 2 values 
(dx 2 /ds = 0). For each galaxy in the sample, we use a red- 
shift given by either spectroscopic or photometric data (see 
Sec." 



2.2 



for details on our sample). 
We define the co-moving luminosity density at the rest- 
frame wavelength A as follows: 

• / \ \ -faint . -mid . 'bright 

Ji(A,Z)—J i -r- Ji + J i — 



$(M z K ,z)fm(M z K ,\)(l + z)dM z K + 

'M 2 =-21.0 

-A/ 2 = -21.0 

§(M z K ,z)m l T l (M z K ,\)(l + z)dM z K + 

'M 3 =-23.0 

rAf 3 = -23.0 

$(Mlc, z)hT t {M z K , A)(l + z)dM z K 

'M 4 = -25.0 

[erg s _1 Mpc _3 Hz -1 ] (3) 

where M^ in Eq. [3] is the rest-frame absolute magnitude 
in K-b&nd at redshift z. The SEDs are given by a function 
Tj(Mjf, A) = L v (with units of erg s _1 Hz _1 ) with i repre- 
senting the different SWIRE SED types. We note that this 
function is dependent on Mjf , since Ti(M^, 2.2) is the lumi- 
nosity per Hz at the effective wavelength 2.2 /im for a galaxy 
with Mk- The fraction of faint, medium, and bright galaxies 
(fi, frii, bi, the three different magnitude ranges) of each of 
the 25 classes is taken into account, according to their M^, 
over redshift. 

The magnitude ranges are denned in rest-frame A'-band 
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Quiescent 



Star-forming 



Starburst 



AGN 



Total non-rejected/rejected 



0.3 235 (29%)/77 (24%) 

0.5 157 (16%)/38 (16%) 

0.7 328 (20%)/59 (13%) 

0.9 144 (14%)/22 (7%) 



554 (69%)/169 (52%) 
756 (77%)/133 (47%) 
1079 (66%)/149 (32%) 
607 (58%)/104 (32%) 



1 (0%)/23 (7%) 

13 (1%)/13 (5%) 

55 (3%)/38 (8%) 

164 (16%)/21 (6%) 



14 (2%)/55 (17%) 

58 (6%)/67 (32%) 

175 (11%)/221 (47%) 

127 (12%)/182 (55%) 



804/324 
984/241 
1637/467 
1042/329 



Table 2. Galaxy SED-type fractions for our galaxy sample after applying the \ d cuts (see Sec. 4.1 1. Numbers are shown for galaxies 
non-rejected and rejected by the cuts, respectively. The total of non-rejected plus rejected galaxies is 5828. This is less than 5986, our 
total number of galaxies, because Le PHARE could not get any fit for 158 galaxies. 
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Figure 7. Galaxy SED-types for the sample after removing the 
worst fits (~ 25% of the total sample, see Sec. |4.1| l. We have 
864 quiescent (in red, 19%), 2996 star- forming galaxies (in blue, 
67%), 233 starbursts (in green, 5%), and 374 AGN galaxies (in 
gray, 8%) from a total of 4467 galaxies. The x-axis describes the 
names of the 25 SED templates from the SWIRE library. 



absolute magnitude using the AB magnitude system accord- 
ing to the different Mi in the integration limits in Eq. [3] The 
faintest magnitude limit corresponds to the faintest galaxy 
in our sample, the medium range is chosen to lead to roughly 
equal numbers of galaxies in each bin in the three magnitude 
bins, and the brightest magnitude corresponds to the bright- 
est galaxy in our sample. This separation is done to take into 
account the fact that for the same fraction of a given SED 
type, the contribution to the luminosity density will be dif- 
ferent depending on luminosity. At any rate, fainter galaxies 
than Mi are too faint to contribute significantly to the EBL 
even if their number density is fairly large. The same is true 
for galaxies brighter than M4: in spite of their high luminos- 
ity their density is not high enough to contribute. 

The function ^(M^,z) in Eq. [3] is the Schechter 
parametrization of the evolving LF as given by C10 in co- 
moving frame. The factor (1+z) comes from the fc-correction 
to account for the change in the definition of the local ab- 
solute magnitude M% with the redshift, i.e., 



k{z) = (1 + z) 



T(M° K ,X) 



(4) 



The co-moving total luminosity density is calculated 
adding the luminosity density from the 25 SWIRE SED 
types, i.e., 



X Quiescent 

X Star-forming 

X Starburst 

AGN 
O Quiescent - SDSS+Galaxy Zoo 
O Star-forming - SDSS+Galaxy Zoo 




0.4 0.6 

redshift 



Figure 8. Galaxy S ED-ty pe fractions from our catalogue (after 
the Xred cu * s : see Sec. 4.1 1 of the different populations versus red- 
shift according to our multiwavclcngth fits. We mark with crosses 
our fractions from z = 0.9 — 0.3. The lines represent the lin- 
ear interpolation that we use to calculate galaxy SED-type frac- 
tions for all redshift: dashcd-red line represents quiescent galaxies, 
dottcd-dashed-blue line represents star-forming galaxies, dotted- 
green line represents starburst galaxies, and solid-gray line repre- 
sents AGN galaxies. The circles at z = 0.1 are fractions computed 
from the SDSS-based sample (see text). We show with a shadow 
area the uncertainties from our lower limit for the errors as well 
as for our v 2 , cut for fits. The uncertainties are around ±0.1. 



jtotal{\z) = ^ jj(\, 



(5) 



We note that the total luminosity density at 2.2 pm 
jtotai(2.2, z) is just the integral of the C10 LF. 

The quantity defined by Eq. [5] gives us an estimate of 
the total amount of light emitted by galaxies per unit volume 
at a given wavelength and redshift. 

The history of the SFR density p in the universe is then 
computed using the following approximation, 

p = 1.74 x 10' 10 {j IR + 3.3J28oo)/£ [M yr^Mpc" 3 ] (6) 

where jm is the total bolometric infrared luminosity density 
integrated from 8-1000 /im, J2800 is the luminosity density 
at 0.28 pm, and L Q - 3.839 x 10 33 erg s" 1 the solar bolo- 
metric luminosity. This equation is taken from |Wuyts et ah] 
(2009), who add the UV and IR contributions (unobscured 
plus obscured), using calibrations for the local universe by 
|Kennicutt|jl998| > and a Salpeter IMF ( |Salpeter|1955) . 

If Eq. [5] is integrated from some redshift 2 to z max — 4 
(up to where the LF is given), the EBL flux seen by an 
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0.3 0.5 



0.7 



0.9 



Total 



faint 


507 


411 


251 


49 


1218 


middle 


255 


452 


899 


530 


2136 


bright 


13 


121 


487 


462 


1113 



Table 3. Number of galaxies (after applying the X^-ed cu ts, see 
Sec. |4.1[ l in every magnitude and redshift bin used to calculate 
the galaxy SED-type fractions in Fig. [9] 



observer at redshift z, due to the radiation emitted from 
Zmax down to z is obtained, 

,2 



A/a(A,z) = 



4ttA 



jtotal[X(l+z)/(l+z'),z' 



(It 



dz 



dz' 

nWnvV- 1 ] (7) 



This is what we call co-moving EBL spectrum, which is 
given in intensity units. The factor dt/dz' takes into account 
the assumed cosmology (e.g., |Peebles||l993| l, and is given 
explicitly by 



tit 



dz' 



1 



H {1 + z')y/Q m {l + z') 3 +n A 



(8) 



with Ho, f2 m and Qa given by the parameters of the ACDM 
cosmology, exactly the same used by C10 for the LF. 

In our approach, it is possible to directly calculate the 
contribution to the EBL from all redshift bins, as well as 
the evolution of the EBL spectrum with redshift and the 
processes related to this evolution, by sources of all the 25 
SED types considered. 



4 RESULTS 

4.1 Galaxy SED-type fractions 

As explained in Sec. [3] the Le PHARE code is used to fit 
every galaxy in our sample to the 25 SWIRE templates. 
For clarity, we will compress in our discussion (but not in 
our calculations, where they will remain independent) the 
25 SED-types in the SWIRE library to four groups: quies- 
cent, star-forming galaxies, starbursts, and AGN galaxies. 
We choose this nomenclature to clarify that our classifica- 
tion is multiwavelength-SED based, and not morphological. 

We note that the fitting procedure is relatively sensi- 
tive to the errors on the photometric measurements leading 
to uncertainties in the galaxy SED-type fractions of ±0.1. 
For our model we set a lower limit of 6% to all the photo- 
metric errors. The effect of different treatments of errors in 
the photometry is discussed thoroughly in Sec. |6.1| and it is 
shown in this section the uncertainties due to this effect on 
our galaxy SED-type fractions and on the EBL estimation. 

To avoid accounting for bad fits, which do not correctly 
describe the galaxy photometric data, a cut in Xred — X 2 1 n 
is applied, with \ 2 given by Le PHARE (Eq. m| and n de- 
grees of freedom (bands with detections). We have checked 
carefully that Xred Sj 30 is a good value for quiescent, star- 
forming and starburst galaxies, but AGN galaxies are sys- 
tematically worse fits, probably due to the fact that there is 
a large range in AGN SED shapes due to multiple emission 




1.2 0.3 0.4 0.5 0.6 0.7 O.i 
redshift 



0.9 1.0 



Figure 9. Galaxy SED-type fractions (after applying the X^ed 
cuts, see Sec. |4.1| of the different populations (the lines are the 
same as in Fig. [Si versus redshift in the three different magnitude 
bins defined in the text for Eq. [3] See Table [3] for number of 
galaxies in every magnitude and redshift bin. 



components which cannot be easily encapsulated in a few 
templates. 

Fig. [6] shows some examples of good and bad fits in 
for the four different main galaxy types with low Xred m 



the top row, some fits around \ r 



10 in the second row, 



other fits around Xr 



30 the third row and some very 



bad fits (with very high Xred) m the bottom row. Due to 
the fact that AGN galaxies are systematically worse fits, 
two different cuts depending on the galaxy-SED type fitted 
are used for our model. These values are Xred ^ 30 for qui- 
escent, star-forming and starburst galaxies, and Xred ^ 10 
for AGNs. As for the uncertainties on the photometric er- 
rors, we show the uncertainties due to these cuts for the 
galaxy SED-type fractions and the EBL and discuss them 
in Sec. [Q 

After applying these cuts there are still 4467 galaxies 
remaining, i.e., ~75% of the original sample. Fig. U\ shows 
a histogram of the galaxy SED-types in the total sample 
after the cuts, and the classification (shown with different 
colours) of the four main galaxy groups considered in this 
discussion. We find 19% quiescent, 67% star-forming galax- 
ies, 5% starbursts, and 8% AGN galaxies. 

A bimodality between quiescent and star-forming galax- 
ies is clearly found. Most of the quiescent galaxies are 
^ 5 Gyr old, late-type elliptical (E115 and E112, according to 
the SWIRE classification). The bulk of the star- forming pop- 
ulation is late-type spirals with the PAH region measured 
using Spitzer data (Spi4, according to the SWIRE classifi- 
cation). In the starburst-like galaxies case, the Arp 220-like 
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Figure 10. Comparison between our estimation of the local luminosity density (black line) and observational data from different surveys: 
12, 25, 60, 100 fira from|Soifer fc Neugebauer| < |l991| l; if -band from|Kochanek et al.||200l) l; u, g, r, i, z, if-ban d from|Bell et al.|l|2003[ ); 
FUV, NUV from |Wyder et a.l.| ( |2005| ); 850 jtm from |Serjeant fc Harrison| ||2005[ l; bj, rf, J, H, A"-band from | Jones et al.| (j2006|) ; 12 , 
25, 60, 100, two different analysis for 170, 800 /im fro m |Takeuchi et al.|j2006| (two different analysis); 8 (an from |Huang et al.| | |2007[ ); 
B-band from|Driver et al,|||2008|l and|Cameron et aT](|2009^; and u, g, r, i, z from |Montero-Dorta fc Prada| )|2009|) , 



galaxies are dominant. The AGN-likc population is clearly 
dominated by Seyfert-type galaxies, especially type II. 

Table [2] and Fig. 15] show the galaxy SED-type fractions 
for four different redshift bins up to 2 = 1, where we have 
chosen bins of Az—0.2 for statistical reasons. The shadow 
regions are the uncertainties due to the lower limits on the 
photometric errors for the catalogue and for the Xred cuts. 
This region is calculated changing the lower limits from 1- 
10% in steps of 1% and applying extreme cases for the cuts 
for every lower limit. The boundaries from these calculations 
lead to the shadow regions. The fractions adopted for the 
model are marked with crosses and wider lines. We observe 
that the fraction of quiescent galaxies increases by a factor 
~ 2 from z ~ 0.9 — 0.3, while the star-forming fraction keeps 
roughly constant for the full redshift range peaking at z — 
0.5. Starburst-type galaxies decrease very quickly from z ~ 
0.9 and reach almost at 2 ~ 0.5. On the other hand, the 
AGN-type fraction is roughly constant from z ~ 0.9 — 0.7, 
and then decreases to 0.02 at 2 ~ 0.3. This result should 
not be considered a complete picture of the evolution of 
the galaxy populations in the universe since these fractions 
depend on the color-magnitude limits of the survey as Fig. [2] 
shows. But what it is certainly described is the population 
of galaxies that contribute the most to the EBL around the 
knee of the LF (the middle and bright region of the LF, see 
Fig. [4}. 

Fig.[9]shows the galaxy SED-type fractions for the three 



different magnitude bins defined in Eq. [3] and explained in 
the previous section, Sec. [3] These are the galaxy SED-type 
fractions used directly in Eq. [3] to calculate the luminosity 
densities. Table|3]lists the number of galaxies in every mag- 
nitude and redshift bin used to estimate the galaxy SED- 
type fractions showed in Fig. [9] 

The galaxy SED-type fractions are extrapolated to 
lower redshift (2 ~ 0) by using results from Goto et al 



(20031, that use data from the Sloan Digital Sky Survey 



(SDSS). They give galaxy fractions according to a morpho- 
logical classification. They are converted to SED classifica- 
tion by using two different observational works, using the 
Galaxy Zoo catalogue from SDSS data, on the abundances 
of blue-elliptical galaxies (/;, e = 5.7 ± 0.4%, Skibba et al. 
|2009[ ) and red spirals (f r3 ~ 25% |Schawinski et al.||2009|). 



Utilizing these works we calculate galaxy SED-type frac- 
tions for the local < 2 < 0.2 universe as follows: in Fig. 12 



and 15 of Goto et al. (20031 we see morphology-density and 



morphology-radius relations respectively. From the bin with 
the largest number density in any of those figures, we have 
the fractions of galaxies with early (~ 14%), intermediate- 
type (~ 26%), early-disc (~ 35%), and late-disc (~ 25%) 
morphology. The fractions of ellipticals are the fractions of 
early galaxies, f e u ~ 14% and the fraction of spirals are the 
intermediate-type, plus the early-disc, plus late-disc galax- 
ies, fspi ~ 86%. From the Galaxy Zoo papers (Skibba et 



al. 2009 



Schawinski et al. 2009 1 these fractions are esti- 
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Figure 11. Comparison between the calculated luminosity densities versus redshift for different spectral bands with observational data 
(solid-black line, for our fiducial extrapolation; dashed-black line for our high- starburst extrapolation for the galaxy SED-type fractions 
for z > 1; see Sec. |4.1[ |. We also show as dot-dashed-orange line the model from [Somerville et al,|2010| Upper-left panel: rest-frame UV 
at 0.28 /im and data from [Gabasch et al.|(2006[l and |Dahlen et a,l.| ( |200T| |. Lower-left panel: rest-frame i^-band and observational data 
from |Arnouts et a.l.| ( |2007p and |Barro et al. | |2009| . It is important to note that this is just the integral of the LF by C10 between M\ and 
M4 in Eq. |3| U pper-right panel: rest-frame B-band and observational data from a co mpilation from|Faber et al. (|2007[l from these works: 
|Norberg etaT] l|2002[| , |Bell et al.| | |2003[l,|Blanton et al.|||2003a| l, [Gabasch et a!1 ( |2004[ l, [E)ahlen et al.| J2005) , arid Ilbert et"aL] | |2005| l. Data 
from |Marchesini et al.|(2007} and |Cameron et al.| | |20 09) arc plotted as well. Lower-right panel: integrated IR from 8-1000 /im data from 
|Rodighiero et al.|j2010l and the phenomenological estimations by Franceschini et al. (2008). 



mated for the local universe according to Eq. [9] and Eq. |10| 
fquies ~ 35% of quiescent and f s f ~ 65% of star-forming 
galaxies. 

Jquies — Jell \Jell X Jbej ~r \Jspi X Jrs) — 

= 0.14 - (0.14 x 0.057) + (0.86 x 0.25) = 0.35 (9) 

fsf = Jspi ~ [fspi X Jrs) + (jell X Jbe) 

= 0.86 - (0.86 x 0.25) + (0.14 x 0.057) = 0.65 (10) 

We have to keep in mind that these numbers are calcu- 
lated from a different sample and a direct comparison with 
our sample may be not accurate. Note as well that our defi- 
nition for quiescent and star-forming is not exactly the same 
as that the red and blue classification from Galaxy Zoo, but 
very similar. Some of our very early-type star-forming galax- 
ies are red according to that classification, but the results do 



not change much because of the fewer number of these galax- 
ies. In the opposite direction to this effect we note as well 
that Le PHARE prefers to fit some early-type star-forming 
galaxies as late-type red galaxies due to their bluer optical 
colours but very little dust re-emission, if any, according to 
the SWIRE templates. 

To be able to compute the local EBL with accuracy, 
as well as its evolution out to the redshifts of the most 
distant objects detected by ground-based VHE 7-ray tele- 
scopes, i.e., z < 0.6 ( Albert et al.||2008 l, we would need to 
extrapolate the galaxy SED-type fractions to 2 > 1. It is 
expected that the local EBL has contributions from these 
larger redshifts, although the behavior is different for the 
optical/near-IR and the far-IR due to the spectral region 
where the different populations contribute. 

For the high-redshift universe (2 > 1, where there are 
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no galaxies in our sample) two different cases are consid- 
ered for the evolution of the galaxy SED-type fraction. It is 
shown that our results are not changed significantly except 
in the far-IR by these two choices. For the redshifts less than 
those of the most distant known 7-ray sources, and redshifts 
where future sources are likely to be found in the near future 
by Imaging Atmospheric Cherenkov Telescopes (IACTs), it 
is found almost no change in the EBL even with a fairly 
large adjustment in the evolution of galaxy-SED type frac- 
tions. This is discussed in Sec. 16. II and here we show the un- 
certainties of the EBL and others quantities calculated due 
to these assumptions. The fiducial choice is to keep con- 
stant the fractions computed for our highest redshift bin. 
This choice is made for simplicity, due to the difficulty in 
the multiwavelength classification of distant galaxies with 
current instruments. But we do note that there is strong 
evidence from several observational re sults by|Reddy et al.| 
( [20051), |Perez-Gonzalez et al.| p008| ), |Taylor et al.| ( |2009f 
and Wuyts et al. ( 2009 1 which suggest no further evolution 
at higher redshifts of the quiescent population. All these 
independent works claim that the fraction of distant non 
star-forming red objects in the high-redshift universe keep 
constant around 24-33% of the total number of galaxies up 
to 2 > 2.5. We find at 2 ~ 0.9 around 14% of quiescent 
galaxies, which it is kept constant for higher redshifts. We 
note here the red-galaxy incompleteness for DEEP2, imply- 
ing that our fractions might underestimate the actual num- 
ber of mainly quiescent galaxies in the faint-end of the LF (as 
seen in the very low number of quiescent fractions in Fig.pjl), 
due to the difficulty for the DEEP2 survey to characterize 
faint-red galaxies for z > 0.8. The impact of this effect is 
decreased by taking into our catalogue galaxies with photo- 
metric redshift. In any case, there are no consequences for 
the EBL results as previously discussed in Sec. |2.2| because 
the bulk of the light comes from the region of the LF around 
Li, where we are basically complete. 

As alternative approach, we choose to increase linearly 
with redshift the starburst-like fraction from our calculated 
16% at 2 = 0.9 up to 60% at z = 2, while decreasing at 
the same rate the quiescent and star-forming galaxies. The 
weight of every one of the 25 SWIRE templates is changed 
in the same proportion. The fractions are kept constant at 
z = 2 for z > 2. This approach is called high-starburst and 
it is used to determine a likely upper limit on the EBL at 



long wavelengths (see Sec. 4.3 1 



4.2 Luminosity densities 

The local galaxy luminosity density is shown in Fig. [To] cal- 
culated using Eq. [5] The solid-black line is computed from 
the sum of the contributions from all the 25 SED-types. 
An excellent agreement is found with the observational data 
from independent surveys over all wavelengths. We note that 
our different assumptions for the high-redshift fractions lead 
necessarily to the same result because this is the light emit- 
ted at z — 0. 

Fig. [Tl] shows the evolution over redshift of the lumi- 
nosity densities at different wavelengths for both of our ex- 
trapolations for the high-redshift fractions. We show with 
dot-dashed-orange line the galaxy formation SAM predic- 
tion by SGPD10 and postpone the discussion to Sec. |6.3| 
The upper- left panel shows the rest-frame 0.28 /im, which is 




1.0 1.5 

redshift 



Figure 12. Comparison between the calculated star formation 
rate density computed using Eq.jfjlfor a Salpeter initial mass func- 
tion, the prediction using the same equation from the luminos- 
ity densities provided by |Somerville et ah] | |2010| | (red-solid line) , 
and some observational data from different estimators shown in 
the legend. The compilation of data points is taken from |Perez-| 
I Gonzalez et al. (2008). The solid and dashed- black lines are from 
the different extrapolations for the galaxy SED-type fractions for 
z > 1 (see Scc.[4~T||. 



in good agreement with the observational results by Gabasch 
|et al.| ( |2006] ) and |Dahlen et al.|p00"7| ) for 2 < 1.5, and some- 
what lower between z = 2—1.5 than the IDahlen et al.| 



(20071 data. This quantity is also directly related with the 
SFR density through Eq.[6] The upper-right panel shows the 
rest-frame _B-band, which is in good agreement with some 



observational results, such as|Norberg et al. ( 2002 1, Gabasch 



et^([2004|),|Ilbert et al.| ( |2005[ ); but around a 15-20% higher 
than|Faber et al. ([2007 1 . At 2 > 1 we are a factor ~ 2 higher 



than the data by |Dahlen et al.| ( |2005[ ). The lower-left panel 
shows the rest-frame K-band, among some observational re- 
sults by |Arnouts et al.| ( [2007) and |Barro et al.| ( |2009"| . The 



lower-right panel shows the evolution of our calculated total 
IR luminosity over redshift, that given by the FRV08 model 
and from observations by Rodighiero et al. (20101. We note 



a general good agreement with these data (we are a a factor 
1.5 higher around 2 ~ 1) for our fiducial extrapolation of 
SED-types beyond 2 = 1, but we are predicting a higher 
luminosity density for the high-starburst assumption. The 
agreement with FRV08 is pretty good, except for the low- 
est redshifts. The total IR luminosity is also directly related 
with the SFR density through Eq. [6] 

4.3 Star formation rate density history 

Fig. 1 12| shows the history of the SFR density of the universe 
computed from our modeling using Eq. [6] It is also plotted 
the prediction using the same equation, from the luminosity 
densities provided by SGPD10, and a compilation of ob- 
servational works from |Perez-Gonzalez et al.| ( |2008[ ) using 
different estimators, assuming a Salpeter stellar IMF. We 
are aware that this IMF is not as good description of the 
observations as other IMFs such as |Chabrier] ( |2003[ ), but 
we are concerned here on showing a comparison with the 
compilation of SFR data, which is given by a Salpeter IMF. 
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Figure 13. The solid-black line is the extragalactic background light calculated by the fiducial extrapolation of the galaxy SED-type 
fractions for z > 1. Empty symbols are direct measurements: 0.3, 0.555, 0.814 /im by |Bernstein| (2007| l ; 1.43, 1.53, 1.63, 1.73, 1.83, 1.93, 
2.03, 2.14, 2.24, 2.34, 2.44, 2.54, 2.88, 2.98, 3.07, 3 .17, 3.28, 3.38, 3.48, 3.5 8, 3.68, 3.78, 3. 88, 3.98 pn by|Matsumoto et al7| < |2005^ using 
IR TS; 1.25, 2.2 jim (slightly shifted for clarity) by |Cambresy et al.||200l| |; 2.2, 2.5 fim by |Gorjian, Wright fc Chary|||2000|j; 60, 100 (im 
by |Finkbeiner, Davis fc Schlcgcl ( 200 0] all these using DIR BE; 65, 90, 140 (slightly shifted for clarity), 160 fim by|Ma,tsuura et aT7|||2010[| 
using AKARI; 100, 14 0, 240 fim by |Lagache et al.| ( |2000| ); 140 (slightly shifted for clarity), 240 /im by |Schlegel, Finkbeiner fc Davis] 
(119981; 140, 240 fim b y|Hauser et al.|||l998[ l all these using FIRAS. Filled symbols are galaxy-count data, usually considered lower limits: 
0.1530, 0.2310 Mm by |Xu et al.| | |2005| l using GALEX; 0.1595, 0.2365 fim by |Gardner, Brown fc Ferguson| 12000] using HST and STIS; 
0.36, 0.45, 0.67, 0.81, 1.1, 1.6 (slightly shifted for cl arity), 2.2 jim (slight ly shifted for clarity) by |Madau fc Pozzetti| j2000| u sing HST 
and ground-based telescopes; 1.25, 1.60, 2.12 /im by |Keenan et al7| | |2010[ | using Subaru; 3.6 /im by |Levenson fc Wr ight (20081; 3.6, 4.5, 
5.8, 8.0 fj,rn by |Fazio et al.| ( |2004fr with a reanalysis of the last point by Franceschini et al. 2008 all these using IRAC; 15 (im b y|Metcalfe| 
|et al.| |2003l l using ISO; 15 /im by |Hopwood et al.| (|2010|) using A KARI; 24 /im by |Papovich et al.|(|2004} an d [Chary et al.| ( |2004[ l; 24 
(sl ightly shifted for cl arity), 70, 160 (im by |Bethermin et al.| | |2010"|) using MIPS; 71.4 /Im by|Frayer et al.||]2006]l using MIPS; 10 0, 160 /im 
by |Berta et al.| l |20To"| l using Herschel. The coloured-solid lines l|Aharonian et al.|2006||Mazin fc Raue|2007||Albert et al.|2008] > are upper 
limits from 7-ray astronomy using different blazars (see Sec . |5| for details) . The dot-dashed-blue line, and the dashed-red line are the 
predictions from the models by Franceschini et al. (2008) and |Gilmore et aL] | |2010] l, respectively. Uncertainties in the our EBL estimation 
are shown with a shadow area. These EBL uncertainties include the uncertainties in Schechter parameters of the LF by Cirasuolo et 



al. 



1 2010), photometric errors in the galaxy catalogue, \ j cuts applied and extrapolations of the galaxy SED-type fractions for 2 > 1 
The envelope of the shadow region within the dashed line at wavelengths above 24 (im shows the region where there is no 



photometry in our galaxy catalogue. The EBL uncertainties are thoroughly discussed in Sec. |6.1| 



The data from z — 3 — 1.5 are roughly reproduced. Our 
results are in agreement within errors with the upper data 
envelope from z — 1.5 — 0.7. We systematically predict a fac- 
tor ~ 1.3 higher SFR than the observational data between 
z — 0.7—0. For the high-starburst assumption a considerably 
higher SFR density is estimated. This high-starburst case is 



motivated by the increasing star formation rate density to 
z ~ 2 in Fig. |12[ and the increasing specific star formation 



rate to z ~ 2 ( |Reddy et al.||2006[ |Daddi et al.||2007[ ). But 
Fig. [12] also indicates that our high-starburst is an extreme 
assumption. 

We want to call attention to the large uncertainties on 
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Figure 14. Extragalactic background light (EBL) in co-moving frame predicted by our model at different redshifts for the two assumptions 
for the extrapolation of the fractions for z > 1 (see Sec. |4.l[ l. The contribution to the EBL from quiescent (red-dashed line), star-forming 
galaxies (bluc-dottcd-dashed line), starbursts (green-dotted line), and AGN galaxies (gray-dotted-long-dashed line) to the fiducial model 
are shown. For comparison, the predictions from other models are shown using magenta-dashed line for Franceschini et al. (2008) and 
orange-dot-dashed line for|Gilmore et aX](|2010). 



Quiescent Star-forming Starburst 



AGN I bol [nW m- 2 sr- 1 ] 



0.0 4.71 (7%) 39.70 (57%) 20.45 (30%) 4.41 (6%) 69.26 

0.2 3.86 (5%) 38.96 (54%) 24.54 (34%) 5.25 (7%) 72.60 

0.6 2.35 (3%) 31.98 (44%) 31.94 (44%) 5.77 (8%) 72.05 

1.0 1.46 (3%) 21.66 (38%) 28.97 (51%) 4.36 (8%) 56.46 

2.0 0.51 (3%) 6.46 (36%) 9.87 (54%) 1.34 (7%) 18.18 

Table 4. Contribution from the different galaxy populations to the bolometric intensity of the extragalactic background light at different 
redshifts in co-moving frame as defined by Eq. 1 1 1 1 to the fiducial extrapolations (see Sec.|4.1[l. 



the observational data estimates for the SFR for all red- 
shifts. These uncertainties are especially important for the 
higher redshifts, mainly because local calibrations are used 
in the estimations, and also the uncertainties of the correc- 
tions due to dust absorption. The same is true for Eq. [5] 
which is calibrated using observed local galaxy properties 
and these might indeed evolve in redshift. 



4.4 Extragalactic background light 

The local EBL (z = 0) estimated using our method is shown 
in Fig. [13] The solid-black line is the EBL calculated by our 
fiducial moderj using Eq. |7J Observational data from di- 
rect measurements (empty symbols) and from galaxy counts 



Intensity files at different redshifts are publicly available at 



http://side.iaa.es/EBL/ 
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Figure 15. Buildup of the extragalactic background light (EBL) 
at different wavelength normalized to z = 0.1. For example, ac- 
cording to the fiducial model ~ 70% of the local EBL at 2.2 fim 
comes from z < 1, but only ~ 40% of the local EBL at 240 /im. 



(filled symbols) are plotted. It is usual to consider data from 
galaxy counts as lower limits. We find a very low background 
from UV to mid-IR, along the lower limits from galaxy 
counts. In the UV our model is lower than the | Gardner^] 
Brown & Ferguson ( 2000 1 data, but we consider these data 



suspect, due to very poor statistics on their number counts 
at the faintest magnitudes and the fact that they are sys- 



tematically higher than the UV data from GALEX (Xu et 
|al.||2005| ), an experiment with higher sensitivity and better 
statistics. 

In the mid-IR region between 7-15 /jm our results are a 
factor ~ 1.2 higher than other models. A lower background 
than FRV08 is estimated from 15-50 /J,m by a factor as large 
as ~ 1.5. Our results are still compatible with the limits from 
galaxy counts. On the contrary, we predict about the same 
far-IR light than FRV08 and a factor ~ 2 — 3 larger than 
GSPD10, higher than the galaxy counts and in very good 
agreement with most of the direct measurements. The high 
flux we predicted in the far-IR (in comparison GSPD10) is 
a characteristic of the SWIRE galaxy SEDs we use, given 
by the GRASIL code which is used to calculate the far-IR 
emission, and the relation between the near-IR and the far- 
IR in the templates. 

In the same figure, we also plot upper limits using solid- 
colour lines from 7-ray attenuation studies. The cyan and 



yellow solid lines by Mazin & Raue ( 2007 1 were computed for 



the so-called realistic and extreme cases, where the authors 
considered different upper limits for the spectral slopes of 
VHE emission from blazars of E~ 1,5 ( Aharonian et al.|2006 



Albert et al 



2008 1 and E 2/3 , respectively. Our calcula- 



tion is compatible with the upper limits from the extreme 
case, but marginally disagrees with the realistic case for the 
largest wavelengths. We will discuss these issues further in 
Sec.El 

Fig.[l3]also shows the uncertainties in our modeling due 
to the uncertainties on the Schechter parameters of the LF 
given by C10, the errors in the photometric catalogue, as ex- 
plained at the beginning of this section, the uncertainties on 
the Xred cu t applied, and uncertainties due to the extrapola- 
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Figure 16. Buildup of the extragalactic background light (EBL) 
at 24 (im obtained from different phenomenological models, nor- 
malized to 2 = 0.1, compared with the Spitzcr/MIPS data by 



|LeFloc'h et al.| p009). For example, according to our fiducial ex- 
trapolation (see Sec. |4.1| , about 75% of the local EBL at 24 (im 
was already in place at z ~1.5. Uncertainties in the modeling 
are shown with a shadow region (see Fig. |13| > - The curve from 
Franccschini et al. (2008) has been calculated by us from their 
published EBL densities. 



tions for the galaxy-SED types for z > 1. All the possibilities 
are calculated and the extreme cases are plotted. The un- 
certainties from the UV up to the mid-IR are dominated by 
the errors in the photometry and the cuts. The directions 
from both effects are different: the uncertainties from the 
photometry are below the fiducial model, the uncertainties 
from the cuts are above it. In the far-IR the uncertainties in 
the extrapolations to Z > 1 dominate with a factor ~ 1.5. 
These effects will be thoroughly discussed in Sec. |6.1| 

The evolution of the EBL is important to account for 
the history of the galaxy emission and the processes in- 
volved, as well as to properly calculate the attenuation for 
VHE 7-rays for the high-redshift universe. We show, in 
Fig. |14| the co-moving intensity level of the EBL for differ- 
ent redshifts, the contribution to the EBL at those redshifts 
from the four main SED groups to our fiducial extrapola- 
tion, and the predictions for the EBL by other models. In 
Table [4] we quantify this evolution, where the bolometric 
intensity is defined according to Eq. |11| i.e., 



hoi 



vl v d\uv 



(11) 



We should note that the starburst population con- 
tributes 54% to the co-moving bolometric EBL at Z — 2, 
but only 30% for the local universe. We note as well that 
the far-IR peak in the SED is higher relative to the near-IR 
peak at these redshifts; this is due to the fact that a large 
fraction of the energy radiated from starburst systems is at 
far-IR wavelengths. We also note that the total bolometric 
intensity peaks at z ~ 0.6 — 0.2, because the far-IR peaks at 
higher energetic wavelengths there as seen in Fig. |14| 

Another important observable is the buildup of the local 
intensities for different wavelengths. This is the fraction of 
the local EBL at a given wavelength that was already in 
place at a given redshift. This is shown in Fig.[l5]for several 
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wavelengths. As an example, we see that ~70% of the local 
EBL at A = 0.445 /jm and 2.2 /im comes from z < 1, 50% of 
the EBL below ~180 /im was already in place at z — 1, but 
it is only ~ 40% at 240 /im. It is significant that the EBL 
at shorter wavelengths mostly come from sources at much 
lower redshifts than the larger ones (sec Lagache, Puget & 
Dole|2005|>. 



Fig. |16| shows a comparison between the EBL buildup 
for our model, FRV08, GSPD10, and the observational work 
by |LeFloc'h et al.||20C)9] ) based on data from MIPS at 24 rim 
up to 2 ~ 1.5 in the COSMOS field. The main contribu- 
tion to the EBL at 24 /urn comes from star-forming and 
starburst-type galaxies. This region of the SEDs is highly 
dependent on the non-smooth PAH features. We observe a 
general agreement, but reaching a factor 40% difference at 
z ~ 1.2 for the fiducial extrapolation. The uncertainties here 



are large (see Sec. 6.1 1 



5 GAMMA-RAY ATTENUATION 

5.1 7-ray attenuation from this EBL model: 
theoretical background 

The EBL has important implications for the interpretation 
of data taken using recent VHE experiments (the Fermi 
satellite, |Gehrels fe Michclson 1999; and IACTs, such as 



MAGIC, VERITA S and HESS; |Lorei^||2u04||Weekes et al.| 
2002 Hinton|2004 respectively), due to the photon-photon 
pair production between 7-ray photons traveling across cos- 
mological distances and EBL photons (sec Nikishov 1962 
Gould fc Schreder|1966 l. 

Blazars are an important source of extragalactic 7-ray 
emission and have become a relevant tool for indirectly mea- 
suring the EBL. These objects are believed to be an extreme 
category of AGNs. Their emission, which occurs at all wave- 
lengths of the electromagnetic spectrum, comes from super- 
massive black holes (with masses Js 10 7 Mq) swallowing 
matter accreted from their surroundings. In general, AGNs 
are characterized by a beamed emission perpendicular to the 
accretion disc known as jets, which are pointing toward us 
in the case of blazars. 

The current theoretical models for the emission by this 
class of objects are of two kinds: leptonic or hadronic. Both 
models predict a spectrum with two peaks, the first one 
localized from radio to X-rays due to synchrotron radia- 
tion from relativistic electrons (leptonic model), or protons 
(hadronic model). However, the second peak has a differ- 
ent nature. While in the leptonic model it is due to inverse 
Compton (IC) scattering of the same population of electrons 
that produce the synchrotron peak ( Bottcher||2007 1, in the 
hadronic model, nuclear photo-disintegration is advocated 
to explain the second peak ( Sikora et al.(2009 |. Both mod- 
els face serious difficulties in explaining intrinsic (i.e., EBL- 
corrected) VHE power law indices harder than 1.5, and fail 
to explain slopes harder than 2/3. The intrinsic spectrum is 
the spectrum that we would observe if there were no effect 
from the EBL. 

The EGRET satellite observed AGNs in the local uni- 
verse (hence not very attenuated), claiming that all of them 

1.5 in the high 
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Figure 17. Upper panel: Optical depth versus observed energy 
of 7-ray photons for sources at different redshifts (from bottom to 
top z = 0.1, 0.3, 0.6 and 1), due to the extragalactic background 
light computed for our model in solid-black line, for Franceschini 
et al. (2008) in dashed-magenta line, and for |Gilmore et alT| (2010) 
in orange-dot-dashed line. Lower panel: Flux attenuation versus 
observed energy of 7-ray photons for fictitious sources at different 
redshifts (from right to left z = 0.1, 0.3, 0.6 and 1). We have cal- 
culated attenuation for the Franceschini et al. (2008) and |Gilmore| 
|et al.|p010| l models using the EBL data provided by the authors. 
The EBL uncertainties in Fig. |13| are propagated to the optical 
depth and flux attenuation. They are shown here with a shadow 
region. 



This result has been confirmed by the Fermi Collaboration 
(within uncertainties), which has published a catalogue of 
AGNs detected by the Fermi Large Area Telescope (LAT) 



all-sky survey during its first year in operation (Abdo et 
|al.||2010a| . From this observational fact, and the theoreti- 
cal issues above, it is usually conservatively consider that 
no AGN could have an intrinsic VHE spectrum fitted by 
a power law with an index harder than 1.5. Some authors 



such as |Katarzynski et al. (20061, Stecker, Baring & Sum- 
merlin|||2007[|,|Bottcher, Dermer fc Finke| J2008|) and|Aharo- 



— n — — — v 

nian et al.| (|20fJ8J) provide some mechanisms within standard 

physics to reach slopes harder than 1.5, but never harder 
than Tint — 2/3. 

The EBL may be constrained using VHE observations 
of extragalactic sources if their intrinsic emitted spectra 
are known. As mentioned in Section [T] 7-ray photons com- 
ing from cosmological distances are attenuated by photon- 
photon pair production by EBL photons. The cross section 
of this reaction depends of the product shown in the left side 
of Eq. [12] 



with Fi 



have spectral indices E 

energy (30 MeV-30 GeV) regime (|Hartman et al.||1999l 



V^eJT 



> 2m e c 



(12) 
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2m„c 



£th 



E(l- cos 6) 



(13) 



where, in the rest-frame at redshift z, E is the energy of the 
7 photon, e is the energy of the EBL photon, and 6 is the 
angle of the interaction, which defines an energy threshold 
E t h for the EBL-photon energy given in Eq. [13] with m e the 
electron mass. 

The cross section peaks at about twice e t h, which pro- 
duces a peak in the interaction at A [/im] = 1.24_B [TeV]. 
From this property, a 7-ray with energy 1 TeV interacts 
mainly with a photon of the EBL with wavelength ~ 1 /xm. 
The details may be found for example in |Madau fc Phinney| 
{T996|. 

For a given observed spectrum of a source at redshift z 
we can find the intrinsic spectrum by assuming a particular 
EBL model and multiplying by the attenuation factor to 
de-absorb the spectrum using Eq. |14| i.e., 



dF 

dE 



dF 

dE 



exp [t(E, z) 



(14) 



where the subscript obs means observed, int is intrinsic, and 
t(E, z) is the optical depth dependent on the observed en- 
ergy E of the 7 photon for a given EBL photon density and 
redshift, 



t(E,z) = 



ft = 



(£)* 



dfi 



/' 



Een{l + z) 2 



de cr 77 (/3')n(e', z')(15) 



(16) 



with dl'/dz' — c\dt'/dz'\ given by Eq. pH fj, = 1 — cos#, <t 77 
the photon-photon pair production cross section, ft' is given 
by Eq.[l6]and n is the proper numb er d ensity per unit energy 
of EBL photons^] We show in Fig. 17 the optical depth and 
attenuation for sources at z = 0.1, 0.3, 0.6 and 1. 

Since the EBL produces an attenuation of the VHE 
spectra, a mere detection of VHE photons (using some con- 
straint on the intrinsic blazar power spectrum) places an up- 
per limit on the EBL density. Some upper limits have been 
derived by different authors, fitting EBL models to the den- 
sity level where the condition Ti nt = 1.5 is satisfied, building 



ad-hoc EBL models. We plotted those limits in Fig. 13 ( Aha 



ronian et al.|2006| |Mazin fc Raue||2007| |Albert et al.|2008 l 
Each of those upper limits comes from the study of different 
blazars with a different measured energy spectrum. Due to 
the peak of the interaction previously mentioned, each of the 
studies constrains different ranges 011 the EBL. Aharonian 
et al.|d2006b used the VHE spectrum of the blazar 1ES 1101- 



232 at z = 0.186 observed from 0.2-3 TeV. scaling the model 
by Primack et al. ([2001 ) multiplying the total EBL intensity 
by a constant to satisfy the Tint = 1.5 condition. |Albert et] 
at] d2008b used the spectrum of 3C 279 at z = 0.536 ob- 



served from 0.08-0.5 TeV, scaling a slightly modified model 



by Kneiske et al. (2002). Mazin & Raue ( 2007 1 used a compi- 
lation of blazars at different redshifts and observed at differ- 
ent energies, and splines from a grid as EBL densities. They 
make two different assumptions about the maximum Tint 



5 Attenuation files are publicly available at 



leading to two different upper limits (called by the authors 
realistic and extreme). 



We see in Fig. 13 that the fiducial EBL model (here- 
after all the results in this section are discussed for this, un- 
less otherwise stated) is below the upper limits at all wave- 
lengths, except at the largest wavelengths, where slightly 
exceeds the limits from the realistic case by Mazin & Raue 



(20071. This fact is discussed in Sec. 5.2.1 and it is ex- 



plained why we do not consider this a major problem. An- 
other limit not plotted comes from the blazar 1ES 0229+200 
at z = 0.140 ( Aharonian et ah]|2007 l. Its study set a lower 
limit in the slope of the local EBL spectrum between 2- 
10 /an, a ^ 1.10 ± 0.25, to satisfy the limit on AGN's spec- 
tra Tint ^1.5. We remark that they set the limit only on the 
slope, not on the intensity level. We have fitted our model 
in that wavelength range, to a power law oc \~ a getting 
a = 1.19 ± 0.07. Our model is thus compatible with this 
constraint. 

It is also possible to set upper limits on the unknown 
redshift of blazars assuming an EBL model and finding 
the redshift by which the EBL-corrected spectrum satisfies 
Tim = 1-5 ( |Prandini et~aL]|2010| [Yang fc Wang||2010[ ). We 
apply that method to the PG 1553+133 spectrum observed 
by MAGIC ( | Albert et aT]|2007[ ), assume an EBL-corrected 
spectrum given by a power law and find an upper limit at 
z ^ 0.85 ± 0.07 in agreement with the lower limit (z ^ 0.4) 



found by Danforth et al. (2010) using absorption features in 



the Ly-a forest of the blazar. 

As shown in Fig. |17| our EBL model implies about the 
same attenuation as other recent models we compare to over 
all the energy range observed by the current generation of 
IACTs. Larger transparency than the observationally-based 
model by FRV08 is found (roughly a factor ~ 2 in flux, 
but still within the uncertainties) for 7-ray photons with 
energies between ~ 6 — 15 TeV for z ~ 0.1, but a factor 
~ 2 in flux less transparent than the GSPD10 theoretical 
approach around ~ 10 TeV. For the large-redshift case, our 
model predicts about the same attenuation as FRV08, but a 
factor ~ 1.5 more transparency than GSPD10 for sub-TeV 
energies. Note that a small difference in the optical depth 
has large effects on the spectra due to the exponential in 
Eq. |14[ e.g., a factor 1.5 in optical depth leads to a factor 
~ 5 in attenuation. 



5.2 Application to extreme known blazars 

We now proceed to test whether the observed spectra of the 
three most constraining AGNs known in the VHE range due 
to their hard spectrum, or to their large redshift, satisfy the 
condition that the intrinsic spectrum corrected by the atten- 
uation derived with our model has Tint ^ 1-5. We consider 
the blazars: Mrk 501 at z = 0.034 detected by the HEGRA 



system of Cherenkov telescopes in 1997 (Aharonian et al. 
|1999[ with a reanalysis by |Aharonian et al.||2001[), FSRQ 
3C 279 at z = 0.536 observed by MAGIC (Albert et al. 



http://side.iaa.es/EBL/ 



2008b , an d the blazar 3C 66 A observed by VERITAS at 
z = 0.444 ( |Acciari et al.|2009| ), all of them seen in a flaring 
state. All these blazars are plotted in Fig. |18| showing in the 
legends that the condition Tint ^ 1-5 is satisfied. 



© 2010 RAS, MNRAS 000,[T}{25 



18 A. Dominguez et al. 



E 
o 



Mrk 501 reanalysis 

z = 0.034 




1 10 20 

Energy [TeV] 



10 

Energy [TeV] 



10" 



10" 



10 



-9 r 



^ io 



■wIts 10 



10 



-10 



li 



-12 



3C279 
z = 0.536 




measured 

intrinsic, r int -3.78ig;J|±0.{ 



10" 



10" 



-10 



E 
o 



5 IO" 11 



6n|Ct) 

Tare 



0.05 0.1 0.4 0.7 

Energy [TeV] 



IO" 12 
0.2 



: ' 1 


■ 


i^^,^^^ 


3C 66A : 

z = 0.444 


I 


I 


+ measured 
: x intrinsic, r 4rtt =1 .42J;g;f |±0.38 


- 


i 





0.4 

Energy [TeV] 



0.6 



Figure 18. Very high energy spectra measured (blue) and EBL-corrected from the attenuation calculated with our EBL model (using 
the fiducial extrapolation for the galaxy SED-type fractions at z > 1, in red) of three extreme blazars: Mrk 501 observed in very high 
state up to energies larger than 20 TeV (upper left panel, |Aharonian et al.|1999| and a reanalysis of the same data (upper right panel, 
|Aharonian et al.|2001[|, 3C 279 a fla t-spectrum radio quasar with the highest redshift (z = 0.536) ever detected for a VHE 7-ray source 
(lower left panel, | Albert et al.|2008| l and 3C 66A a BL Lac with probably (because its redshift, z = 0.444 is not very secure) the highest 
redshift ever detected for an object of this class (lower right panel, [Acciari et al.||2009[ |. Uncertainties from the EBL modeling as well 
as statistical and systematic errors are shown with a shadow region. The straight red line is the best-fitting power law for every blazar 
with index f\ n t. The first uncertainties in the index are due to the EBL modeling as shown in Fig. |17| and the second uncertainties are 
statistical plus systematic errors in all blazars, except 3C 66A where only statistical errors are shown. 



5.2.1 Mrk 501 

The highest energy bins in this measurement, where it 
is observed a significant deviation from a power law (see 



Fig. 18 upper- left panel), are affected by the far-IR EBL at 



A > 60 fira. This is the region of the EBL spectrum where 
it was found a disagreement with the realistic (but not ex- 
treme) upper limits of Mazin & Raue (20071. The problem 



comes from the very low statistics and high systematic un- 
certainties at such high energies ( jAharonian et al.|1999 l. A 
later reanalysis of the same observation done in I Aharonian 
et al. (20011 accounts for larger systematic uncertainties as 



shown in Fig. [T8| upper right panel. 

This exponential behaviour for the highest energy 
bin was already observed from the first EBL models 



(e.g., Malkan fc Stecker|1998 Primack et al.||1999 Kneiske 



et al. 2002), whose EBL levels were higher than the more 
recent ones. This fact was discussed thoroughly in |Dwek fc| 
Krennrich (20051, and even some exotic explanations such 



as Lorentz invariance violation (Stecker & Glashow 20011 



were proposed. More recent EBL models with a more trans- 
parent universe (such as our model, FRV08 and GSPD10) 
relax such predictions. The solutions to exponential spectra 
and photon pileup could involve widespread problems with 
the photon statistics and systematic uncertainties in the ob- 
servations (as suggests the results from the later reanalysis), 
or new mechanisms extending the normal SSC model, using 
external regions close to the 7-ray source with target pho- 
tons. The EBL uncertainties in the far-IR leading to the 
attenuation uncertainties at these high energies as shown in 
Fig. |18| might contribute to the solution as well. 

Another observed flare with better statistics with the 
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current generation of IACTs up to such high energies as 
~ 20 TeV would be very helpful in constraining these pos- 
sibilities. 



5.2.2 3C 279 

Fig.[l8]shows in the lower-left panel the EBL-corrected VHE 
spectrum for this source. An external photon field providing 
target photons for IC (such as that provided by a broad- 
line emission region) might be necessary to explain the flat 
behavior at the largest energy bins, as discussed in Albert 
et al. (20081. Instrumental systematic uncertainties might 



explain this behavior as well. We note here that our model 
is already matching the lower limits from galaxy counts at 
the wavelengths where 7-ray attenuation with the observed 
energies occurs, and a much lower EBL density than the 
one calculated in this work does not seem realistic. The at- 
tenuation uncertainties from the EBL modeling are too low 
at this redshift and these energies to explain that spectral 
behaviour. 



5.2.3 3C 66 A 

Fig. [18] shows in the lower-right panel the EBL-corrected 
VHE spectrum for this source, whose EBL-corrected slope 
is well within the la limit of the Ti nt ^ 1.5 according to the 
calculated intrinsic index. It is important to note that the 
redshift considered for this object is calculated using just 
one emission line and is thus not very secure (see discussion 
Bramel et al.|2005l, and its attenuation might be indeed 



overestimated if the redshift is lower than assumed. 



5.3 Conclusions on the limits from blazars 

It is concluded from the study of these extreme blazars that 
our EBL is generally compatible with the hardness of the 
EBL-corrected slopes expected from theoretical arguments. 
However, it is clear that a simple SSC model cannot explain 
any flatness at the highest energies of the EBL-corrected 
spectra of either Mrk 501 or 3C 279, which suggests that 
some extension to the model may be necessary such as an 
external photon region, a better understanding of the IACT 
systematic uncertainties, or even a revision of the propaga- 
tion mechanisms mainly through the intergalactic medium 
(see Sanchez-Conde et al.|2009l. 



5.4 Propagation of the EBL uncertainties to the 
7-ray attenuation 

As shown in Fig. [17] the uncertainties in the attenuation are 
dependent on the observed 7-ray energy as well as the red- 
shift: the higher the energy or the redshift, the higher the 
uncertainties in the attenuation. The attenuation uncertain- 
ties shown were calculated from the uncertainties in our EBL 
modeling, which were shown in Fig. |13| explained in Sec. H] 
and will be thoroughly discussed in Sec. |6.1| For sub- TeV 
energies up to around 1 TeV the uncertainties in the flux 
attenuation are never higher than a factor ~ 2 and gener- 
ally lower. The uncertainties in the EBL-corrected spectra 
in this case are dominated by other effects (see indices in 



are around a factor of several. The uncertainties in the EBL- 
corrected spectra up to such high energies due to the EBL 
modeling are considerable. These high uncertainties are de- 
rived from the EBL in the far-IR region due to the very fast 



increment of the EBL photon density (n, see Eq. 15 1 with 
longer wavelengths. Observations of sources at low redshift 
but energies larger than ~ 10 TeV will set constraints on 
these uncertainties. 



6 DISCUSSIONS AND COMPARISON WITH 
SEMI-ANALYTIC MODELS 

6.1 Discussion on EBL uncertainties 

As explained in Sec. [3] we adopt a lower limit to the pho- 
tometric errors higher than those in the AEGIS catalogue. 
Different lower limits are set from 1-10% of the photometric 
measurements. That is, if the error in any band is lower than 
our limit then we set it to the limit. The results are sensi- 
tive to the limit choice in R and I (where the errors in the 
catalogue are the lowest), but not for the other bands. The 
galaxy SED-type fractions change for lower limits 1-6%, but 
there is little change if the level is set higher than 6%. The 
change is mostly for the quiescent and star-forming galaxy 
fractions. If we use the errors in the catalogue without any 
change, we find 10% more quiescent galaxies at z = 0.3 than 
for a lower limit of 6%, which decreases to ~ 3% more qui- 
escent galaxies at 2 = 0.9 than for a lower limit of 6%, as 
shown in Fig. [8] The change is mostly in E112-type galax- 
ies, according to the SWIRE classification. We have inves- 
tigated those quiescent galaxies that change their best fit 
to star-forming galaxies upon raising the lower limit on the 
errors, and find that they are often fitted much better by a 
star-forming SED. In many cases they even have detection 
in MIPS 24, clearly indicating ongoing star formation. On 
the other hand, based on the comparison of our photomet- 
ric measurements to those of other catalogues we estimate 
that any error in the photometry lower than ~ 5% should 
not be considered very reliable. For those reasons we set the 
lower limit at 6% for model. The uncertainties due to this 
are below the fiducial model in Fig. |13| and for the reasons 
stated an EBL in this region should not be considered very 
likely (and therefore, neither is their derived attenuation in 



Fig. 171 



Anot her source of uncertainty accounted for in Fig. [8] 
is the Xred cu t that separates good and bad 



13 



and Fig. 

SED fits. The main change occurs for AGNs, where for a 
relaxation in the cut (from Xred — 10 to 20), the fraction 
can increase by as much as 10% at Z = 0.9, and by a smaller 
fraction at z — 0.3. These changes affect the EBL in the fol- 
lowing way: higher AGN fractions increase the UV as well 
as the mid-IR, while higher quiescent fractions decrease the 
flux at those wavelengths. This effect affects the uncertain- 
ties above the fiducial model and an EBL intensity in this 
shadow region is considered more reliable than the region 
below the fi duci al model (the same for their derived attenu- 



ation in Fig. 17 1. While Xred cu ts do not have an appreciable 



Fig. 18 1. For energies larger than 10 TeV the uncertainties 



effect on the far-IR flux, there is a substantial change arising 
from the choice of extrapolation in SED-types above z ~ 1, 
as we find in the high-starburst assumption. Fig. |13| also ac- 
counts for the uncertainties in the Schechter parameters of 
the LF given by C10 but these are small. 
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Two major potential problems for our modeling might 
be a colour-dependent selection effect and the extrapolation 
of the galaxy SED-type fractions for z > 1. It was already 
shown in Sec. 12.21 that the colour-selection effect is rather 
small. From the fact that most of the light in the EBL comes 
from the knee of the LF around L*, where our sample does 
not suffer any color-dependent selection effect, we do not 
consider this to be a significant problem for our EBL calcu- 
lation. Our estimated galaxy SED-type fractions appear to 
be consistent with works by others as well. For example, our 



results agree with Blanton (20061 and Faber et al. (20071, 



who find roughly no evolution for late-type (blue) galaxies 
from z ~ 1 — within a 10% range, and an increment of the 
early-type (red) population in the same redshift range by at 
least a factor 2. We also highlight that the galaxy SED-type 
fractions that we calculated for the local universe smoothly 
link with our independently-derived results at z ~ 0.3. 

Regarding the galaxy SED-type fraction extrapolations, 
we have considered two rather different approaches which 
basically lead to the same evolving luminosity densities and 
EBL for the optical/near-IR range where 7-ray attenuation 
occurs, as shown in Fig. [TTJ and Fig. |14| respectively. This 
fact is due to the shape of the stellar emission, because the 
contribution to the optical/near-IR peak is very similar for 
quiescent, star-forming and starburst galaxies for a given 
Mk ■ We recall here that the normalization to our model is 
fixed by the rest-frame if-band LF by C10. The only differ- 
ence between our extrapolations is at the far-IR peak, where 
our results are considered to be less robust for this reason, 
as well as for the reasons stated below regarding the SWIRE 
templates and the lack of photometric data. Deeper obser- 
vations by future galaxy surveys will help in characterizing 
the galaxy SED-fractions up to higher redshifts. 

It was also checked how the fractions change if the de- 
tection limit is relaxed from 5a to 3a for the bands where 
there are observations, but no detections. Many more quies- 
cent galaxies than in the 5a case were found, even a factor 
larger than 2, due to the low detection limit on the MIPS 24 
instrument, but generally they are not good fits. 

In our work we have not differentiated between the 
spectroscopic and secure photometric redshifts. This is an 
approximation, and it is necessary to check that this is con- 
sistent with our results. We find that the galaxy SED-type 
fractions derived from both sub-catalogues are clearly com- 
patible and the trends are the same. 

Galaxies fitted to a starburst SED may instead be very 
late- type star- forming galaxies (or viceversa), because both 
SED templates are rather similar in the regions where we 
have data. This may be called mis-typing and its effect is 
expected to be larger for faint galaxies, because the ma- 
jor fraction of faint galaxies are star-forming or starburst 
and not massive quiescent or AGN galaxies. Such small 
galaxies are probably rather metal poor and thus lacking 
dust. Hence, their SEDs are probably more like star-forming 
galaxies rather than starbursts (in agreement with our re- 
sults in Fig. |9l. It is a source of uncertainty in the mid and 
far-IR (underestimating or overestimating light) and might 
explain the excess found in Fig. [TTJ and [12] compared with 
the data. Far-IR data would help in resolving this issue, but 
the number of galaxies with detection in MIPS 70 is rather 
low to make statistical estimations. Herschel data will be 



very useful thanks to its good spatial resolution and deep 
photometry in the far-IR. 

Another source of uncertainties in our model that quan- 
titatively we have not accounted for arises from the use of lo- 
cal SED templates to fit galaxies at z > 0.3. This comes from 
the fact that the SWIRE templates are based on observa- 
tions of local galaxies, and we expect that they become worse 
fits when the redshift is increasing. This problem will be ad- 
dressed by new data from WFC3 on Hubble Space Telescope 
and the next generation of ground and space optical/near-IR 
telescopes such as the James Webb Space Telescope. 

The lowest EBL flux in Fig. [13] is given by the case 
with the highest number of quiescent galaxies and lowest 
number of AGN galaxies, which corresponds to the case of 
using the low errors in the catalogue and our Xred cuts. 
The highest EBL flux occurs with fewer quiescent galaxies 
and the highest fraction of AGN galaxies. This is the case 
with the 6% lower limit for the errors in the photometry and 
without a Xred cut. Using the lower limits from galaxy counts 
in the UV and in the optical we may rule out at > 2a the 
mixing of galaxy SED-type fractions predicting the highest 
fractions of quiescent galaxies and fewest AGNs in Fig. IS] 
We do not consider the VHE observations to exclude the 
models with higher far-IR, because the discrepancy is for 
wavelengths longer than 60 /^m where those limits may not 
be reliable for the reasons stated in Sec. [5] Further VHE 
observations might indeed constrain our galaxy SED-type 
fractions. 

Thus to recap, the EBL uncertainties from the UV up 
to the mid-IR are low enough to recover the spectra of 7- 
ray sources with energies lower than ~ 10 TeV, but the 
EBL uncertainties have to be reduced in the far-IR (for 
neglecting uncertainties due to the EBL modeling) to cor- 
rect higher energy sources. Additional photometry is needed 
there to clearly distinguish between star-forming and star- 
burst galaxies, therefore to reduce the mis-typing, as well as 
a better understanding on the far-IR region of the galaxy 
SEDs at z > 0.3. Characterizing the galaxy SED-type frac- 
tions at 2 > 1 will reduce these uncertainties in the far-IR 
as well. 



6.2 Discussion of the results 

The local luminosity density from galaxies is observationally 
well constrained over all wavelengths from 0.1-1000 /im. As 
shown in Fig. |10| the prediction of the local luminosity den- 
sities is in very good agreement with observational results. 
Fig. [ll] showed the evolving rest-frame UV luminosity 
density as well. These results agree well with the observa- 



tional data by Gabasch et al. (12006 ) and Dahlen et al. ( 2007 1 
within uncertainties, but they are a factor ~ 1.6 below the 
data for z = 1.7 and 2.2. It was also compared in Fig. [TTJ 
the evolution of the rest-frame luminosity in the K-ha,n& 
from our calculations to independent observational works. 
Some disagreement was found that in the case of Barro et 
al. ( 2009 1 might be explained by the fact that they do not 



correct their sample for incompleteness and only consider 
the brightest sources, unlike the LF by C10. Therefore their 
results should be considered as lower limits. The direct com- 
parison in Fig. [TTJ with the rest-frame luminosity density in 
B-band showed that our luminosity in that band is not in 
contradiction with other independent works. We are doing 
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really well reproducing the data from Norberg et al. (20021, 



Gabasch et al. (20041 and Ilbert et al. (20051. We might 



be indeed overestimating the light in this band 15-20% for 



z < 1 according to the data by Faber et al. (20071 and a 
factor ~ 2 for z > 1 according to the data by |Dahlen et| 
|al.| ( |2005[ l, but this latter does not significantly affect our 
results for the local EBL because as we already showed in 



Fig. 15 most of the optical/near-IR light comes from z < 1. 



The comparison with the bolometric IR luminosity density 
with the observational works by FRV08 and Rodigliicro et 
al. (20101 is very good, even though we are a factor ~ 2 
higher around z ~ 1. 

A good agreement was found with the upper envelope of 
the data cloud on the calculations for the SFR history from 
z ~ 1.5 down to the local universe using Eq. [6] (see Fig.|12[), 
using our fiducial extrapolation. According to Magnclli et 
al. ( 2009 1 at around z ~ 1 the main contributor to the star 



formation is the obscured IR contribution, instead of the 
UV. We may be overpredicting some of this obscured IR 
light around z ~ 1 due to the lack of far-IR photometric 
data in our galaxy catalogue that allow a clear classifica- 
tion between late star-forming and starburst galaxies as dis- 
cussed in Sec. |6.1| With the high- starburst extrapolation of 
the galaxy-SED-type fractions was checked that increasing 
the starburst-like population a factor ~ 3 from z ~ 1 — 2, 
we may get a flatter SFR density history up to z ~ 2, but 
even higher than the observational data. This high- starburst 
assumption does not change our general picture of the local 
EBL, but increases the far-IR peak a factor ~ 1.5 (as it was 
considered in Fig.|13|and was shown explicitly in Fig. |14| for 
some other redshifts). 

The EBL calculated in this work is matching the data 
from galaxy counts from the UV up to the mid-IR (see 
Fig. 13 1, except the data found in Levenson & Wright 2008 



Higher intensities than the data from galaxy counts were 
calculated in the far-IR but in agreement with direct de- 
tections. The EBL evolution shown in Fig. [14] is in good 
agreement with FRV08 up to z = 1. At higher redshift our 
results are different in the UV and optical/near-IR. This 
may be due to the fact that FRV08 extrapolate the galaxy 
evolution, while in our model this evolution is entirely based 
on the observed LF by C10 up to z = 4. See Sec. |6.3| for a 
comparison with the results by GSPD10. 

There are some works in the literature where the con- 
tribution from AGN galaxies to the EBL is studied. Accord- 



ing to recent works that focus in the mid-IR (e.g., Silva, 
Maiolino fc Granato|2004| |Matute et al.|2006[ ) this contribu- 



tion should not be larger than 10-20%. This is in agreement 
with our results: we find that the AGN-galaxy contribution 
to the bolometric EBL is 6% for the fiducial extrapolation 
(Table [4} and 13% for the case with the largest AGN frac- 
tion in Fig. 15] For the wavelength range between 1-20 fim 
the AGN contribution from our model is also between 8- 
16%. We estimate that this contribution to the co-moving 
bolometric EBL slightly increases with redshift. 

The EBL buildup was studied in Fig.[l5]and[T6] It was 
found that most of the local UV/optical/near-IR EBL was 
built up at z < 1, while the far-IR EBL was mostly built 
up at z > 1. This result for the far-IR light agrees well with 
the observational work by Devlin et al.| (|2009[) , but disagrees 



in the buildup of the local EBL at 24 /im. These differences 
are due to the fact that a very small change in the mid-IR 
region leads to a very strong difference in this buildup plot, 
and to the mid-IR peak that we get at larger redshifts (see 
Fig. 14 1 due to the shape of the galaxy SEDs. We point out 



with Chary & Pope (20101. In any case, our uncertainties in 



here that the EBL buildup is on how the light is being built 
up, and not about the absolute intensity value. 

We already discussed our results on 7-ray attenuation 
in Sec. [5] 



6.3 Comparison with SAMs 

In this section we compare our EBL estimation against the 
EBL model described in SGPD10 and GSPD10, which is 
based on SAMs of galaxy formation. The comparison for 7- 
ray attenuation has been thoroughly discussed in Sec. [5] We 
notice that slightly different cosmological parameters were 
used for our model and that by SGPD10. The latter uses 
the latest values from WMAP5, which slightly affect all the 
results in the local universe as well as their evolution. 

We already saw in Fig. [TT] the comparison between our 
observational luminosity densities and the theoretical pre- 
diction by SGPD10 for the co-moving luminosity density 
versus redshift in the UV, in the near-IR (K-band), in the 
optical (_B-band) and for the bolometric IR luminosity. We 
note that our K-b&nd luminosity density evolution is given 
exclusively by the C10 LF, because at that band our choice 
of galaxy SED fractions does not affect our results. This 
quantity is above the prediction by SAMs by a factor around 
20% from z ~ 2 down to the present universe. The UV from 
this SAM is above our results for all redshift, except at z ~ 1. 
At z ~ 2 is a factor 4 higher. For the B-band luminosity den- 
sity the agreement is excellent from z ~ 1 down to the local 
universe. For z > 1 SAMs predict a factor of several more 
light than our observationally-based approach. We may see 
the consequences of this for the EBL evolution in Fig. [14] 
for high redshifts where the excess of light has not been di- 
luted by the expansion of the universe. For the bolometric 
IR luminosity SAMs seem to systematically predict at least 
a factor ~ 2 less light than our calculations and the one by 
FRV08. This difference is maximized around 2 ~ 1 up to a 
factor ~ 4. 

Fig. [12] showed a comparison between our SFR density 
estimation and that predicted by SGPD10 as calculated by 
using the Eq. [6] From z ~ 3 — 1 our SFR densities have 
different a behaviour: for our observational model increases 
up to z = 1 and for SGPD10 keeps constant down toz~ 1.7. 
For lower redshifts both models decline down to the local 
universe. 

In general a very good agreement between the local EBL 
from the UV up to the mid-IR predicted by our method and 
the SAM of SGPD10-GSPD10 was seen in Fig. [l3j A factor 
~ 1.5 higher intensity is found in the local UV from SAMs, 
and around the same factor lower intensity around 15 fim. 
For the far-IR peak, the difference comes from the different 
templates used for the dust component in the far-IR, which 
is given by the GRASIL code in the case of the SWIRE 
templates (which we use), and by a interpolation between 
the observed 70 and 160 /im by MIPS in the case of the 
templates used by these SAMs ( |Rieke et al.|2009[ l. 

The agreement on the evolution is very good as well 



the far-IR are very large. Differences up to 40% were found as seen in Fig. 14 even though at high redshift (z > 2) 



© 2010 RAS, MNRAS 000,[T}{25 



22 A. Dominguez et al. 



GSPD10 predicts a factor of several more light in the UV. 
This is due to galaxies within the faint end of the theoretical 
LF at z > 2. We recall that our observational model seems 
to already be overproducing light in the B-band for z > 1.5 



according to data in Dahlen et al. (20051 (see Fig. Ill 



6.4 Overview on the cosmological picture 

The evolutionary path that we have in mind to interpret the 
evolution of the galaxy SED-type fractions is the following: 
AGNs are mostly formed by mergers between galaxies dur- 
ing hierarchical growth of dark matter halos (e.g., Hopkins 
|et al.|[2~009[ ). When the merging galaxies are gas rich (also 
known as wet mergers), inflows of gas are produced leading 
to starburst galaxies and to the mass growth of the central 
black hole. The central black hole activity begins to expel 
the gas. Eventually, the gas is exhausted, switching off the 
AGN. The galaxy continues forming stars as a star-forming 
galaxy until the gas is fully depleted, then becomes a quies- 



cent galaxy ( Hopkins et al.|2008a Hopkins et al.|2008b |. 

It is now a well known observational fact that galax- 
ies are bimodal in some properties such as colours (Strat- 



eva et al. 2001 Blanton et al. 2003b I. They group in two 



different regions in colour-magnitude diagrams defining the 
red sequence and the blue cloud. Galaxies forming stars are 
in the blue cloud. Some galaxies have their star formation 
quenched when they become satellite galaxies in a larger 
halo, they cease to accrete gas, and they join the red se- 
quence. Central galaxies form in the blue cloud, but they 
join the red sequence when they form a supermassive black 
hole and/or their halo mass exceeds approximately 10 12 Mq 
and/or they become satellite galaxies in a cluster. The most 
massive red galaxies cannot have simply be quenched cen- 
tral blue galaxies, since the latter are not massive enough; 
thus they must have been created by mergers without much 
star formation (also known as dry mergers). This effect is 
shown in Fig. Ill taking into account that the Tf-band ab- 
solute magnitude Mk is a good tracer of the galaxy stellar 



mass, as shown in IBrinchmann & Ellis ( 2000 \ 



Massive galaxies today (very bright Mk) form their 
stars first, which is known as downsizing (|Cowie et al. 1996 



Perez- Gonzalez et al. 20081. This initially seemed at odds 



with the hierarchical nature of the ACDM paradigm, in 
which small halos form first and agglomerate into larger 
ones. But the idea that star formation is efficient only in 
dark matter halos with a narrow range of masses naturally 
explains how the phenomenon of downsizing arises: halos 
that are massive today passed through the star-forming mass 
band between 10 8 -10 12 Mo earlier and thus formed their 



stars earlier than halos that are less massive today (Cro- 
ton|2009| [Conroy fc Wechsler|2009|> 

A careful examination of Fig.HJreveals some interesting 
trends. We find that all the oldest galaxies (~ 13 Gyr old, 
E1113) are in the red sequence. However, the younger quies- 
cent galaxies (~ 5 Gyr old, E115) can be found in the red se- 
quence as well as in the green valley (the region between the 
red sequence and the blue cloud) . For the youngest quiescent 
galaxies (~ 2 Gyr old, E112) we find that for z > 0.6 they 
populate the green valley, while for z < 0.6 they belong to 
the blue cloud. All the early-type star-forming galaxies (SO, 
Sa) are in the red sequence. Later-type star-forming galaxies 
such as Sb and Sc start to populate the green valley as well 



as the red sequence. Most of the very late-type star-forming 
galaxies (Sd, Sdm, Spi4) populate the blue cloud. Starburst 
galaxies are mainly in the green valley, but some of them 
are in the bluer region of the red sequence and in the red- 
der region of the blue cloud. The same happens to AGNs, 
but they tend to be in the blue cloud more than in the red 
sequence. 

We note that the increasing rate of quiescent galaxies 
as z declines is roughly the same as the decreasing rate of 
starburst-type galaxies from z ~ 0.9 — 0.7. One possible ex- 
planation would be the direct transformation of starbursts 
(either merger or huge-cold-gas reservoir triggered) directly 
to quiescent galaxies, without an intervening stage of signifi- 
cant star formation. Another explanation is that the charac- 
teristic time in which starburst-like galaxies consume their 
cold gas is the same that in which star-forming galaxies con- 
sume their lower cold-gas reservoir. Thus the specific SFRs 
of these populations are different, but the rate at which star- 
bursts enter the star-forming sequence is the same as the 
rate at which star-forming galaxies become quiescent. From 
z ~ 0.7 — 0.3 the fraction of starbursts is very low, so the 
constant increase of the red sequence is modeled as due to 
AGNs preventing gas from cooling and forming stars. 



7 SUMMARY AND CONCLUSIONS 

A novel, robust and powerful method based on observations 
to derive the evolving spectrum of the extragalactic back- 
ground light (EBL) between 0.1-1000 fim was presented. 
This model is based on the observed rest-frame K-b&nd 



galaxy luminosity function (LF) over redshift by Cirasuolo 
et al.|p010[ ) (C10), combined with an estimation of galaxy 
SED-type fractions based on a multiwavelength sample of 
~ 6000 galaxies from AEGIS. This model has the following 
main advantages over other existing EBL models: transpar- 
ent methodology, reproducibility, and utilizing direct galaxy 
data. The best available data sets are used (ClO's LF and 
the AEGIS galaxy catalogue) observed over a wide redshift 
range. The galaxy evolution is directly observed in the rest- 
frame K-band up to z = 4. Observed galaxies up to z = 1 
from the UV up to 24 /jm with spectral energy distributions 
(SEDs) of 25 different types (from quiescent to rapidly star- 
forming galaxies, and including AGN galaxies) are taken 
into account in the same observational framework. A study 
of the uncertainties to the model directly from the data (such 
as uncertainties in the Schechter parameters of the C10 LF 
and the errors in the photometric catalogue) was done, and 
their propagated uncertainties to the 7-ray attenuation were 
studied. 

A brief comparison with results from other recent EBL 
models is made here: Stecker et al. (2006) estimate a lo- 
cal EBL in the UV a factor of several higher than us, and 
in contradiction with recent 7-ray observations (|Abdo et 



al. 2010b I. A comparison with the FRV08's results was 



throughly presented through Sec. HI In general, our results 
are in good agreement, despite the fact that our modelings 
are different. Finke et al. (2010) have five different mod- 
els based on different parametrizations of the SFR density 
of the universe and IMFs. The local EBL from the UV to 
the near-IR are similar to ours for their models C and E. 



Kneiske & Dole ( 2010 1 claim to model a strict lower limit for 
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the EBL. However our results for the local EBL in the UV 



are lower than the calculation by Kneiske & Dole ( 2010 1, but 



are in agreement with the robust lower limits from galaxy 



2010) is not compatible with the lower limits by Keenan 



counts in the UV by Madau & Pozzetti ( 2000 1 and Xu et 



al. (20051. In the near-IR, the model by Kneiske & Dole 



et al. (20101. Our observationally-based approach was also 



throughly compared with the results from the semi-analytic 
model (SAM) of galaxy formation by SGPD10 and GSPD10 
in Sec. |6.3| Our EBL results are in general in good agree- 
ment at least for z < 2, even though this SAM predicts more 
light (by up to a factor of several) than our observational 
approach in the UV, and a factor ~ 2 — 3 less light in the 
far-IR. 

Our methodology provides a tool for calculating the 
EBL more accurately at the longest wavelengths when a 
better understanding on the far-IR galaxy SEDs, new pho- 
tometry, and deeper LFs at those wavelengths are available 
from the Herschel Space Observatory. 

Two extrapolations of the galaxy SED-type fractions 
to z > 1 were considered, showing that these assumptions 
only affect the far-IR. It was calculated that the popula- 
tion with SED features of quiescent local galaxies increases 
a factor ~ 2 since z ~ 1. The star- forming population re- 
mains roughly constant, while the starburst-like population 
decreases very quickly from around ~ 20% at z ~ 1. The 
AGN-like population decreases slower than the starburst- 
like population from almost 20% at z ~ 1 to just around 2% 
at z ~ 0.3. Data from the future James Webb Space Tele- 
scope will help to determine the galaxy SED-type fractions 
at z > 1. 

A low intensity local EBL (z — 0) was found, match- 
ing the lower limits from galaxy counts up to ~ 30 /an. 
For longer wavelengths, our model predicts higher intensi- 
ties than the data from galaxy counts, in agreement with 
direct measurements. 

Our results are also compatible with all the upper limits 
from 7-ray astronomy according to the standard framework 
for the propagation of VHE photons through the universe, 
even though to account for the highest energies detected by 
Aharonian et al. J1999) for Mrk 501 we have to assume a 



^int < 1.5, appeal to statistical and systematic uncertainties 
on this VHE spectrum, or attenuation uncertainties due to 
uncertainties in the EBL for such high energies as discussed 
in SecJU 

The EBL uncertainties in far-IR leading to attenuation 
uncertainties of a factor of several for energies larger than 
~ 10 TeV needs to be addressed by the current and next 
generation of IR telescopes providing new photometric data 
and a better understanding of the galaxy IR emission. 7-ray 
astronomy may constrain these uncertainties from a better 
understanding of the emission mechanisms at those high en- 
ergies (helped with simultaneous multiwavelength observa- 
tions) and of the instrumental systematic uncertainties. Ob- 
servations aimed to measure photons with energies higher 
than ~ 10 TeV at 2 < 0.3 are encouraged. 

It is worth mentioning that high energy (30 MeV- 
30 GeV) 7-rays are detected by Fermi for z ^ 2.5 from 
AGN s ( |Abdo et aLl|2010a[ ) and for z < 4.5 ( |Abdo et al" 
2009) from GRBs. The reasons for these high redshift de- 
tections include a larger 7-ray flux at lower energies and a 
lower density of EBL target photons that can interact with 



these 7-rays. Understanding the evolution of the EBL at 
UV wavelengths is essential to interpreting observations of 
these high-redshift sources. New observations of AGNs as 
well as the first GRB detection in the VHE range would 
help to make new and stronger constraints on the EBL; see 
|Gilmore, Prada fc Primack| polo) . 

The universe, according to our observationally-based 
model, is more transparent than the estimation from FRV08 
(a factor ~ 2 in flux) for VHE photons coming from low- 
redshift sources (2 ~ 0.1) for energies between ~ 6 — 15 TeV, 
but still the uncertainties here from the EBL modeling are 
large (a factor of several). The same attenuation than FRV08 
is estimated for other energies. For VHE photons coming 
from larger-redshift sources (2 ~ 1), roughly the same at- 
tenuation as FRV08 is estimated. Here the attenuation un- 
certainties (for energies available to 7-ray telescopes) due to 
the uncertainties on the EBL modeling are low in compar- 
ison with other effects. At these redshifts the uncertainties 
on the EBL-corrected spectra are dominated by instrumen- 
tal systematic uncertainties. We may conclude that it is not 
expected to observe any such high redshift (2 ~ 1) multi- 
TeV 7-ray photons from blazars with the current or even 
next telescope generation such as the Cherenkov Telescope 
Array (CTA,iDoro 2009) or the Advanced Gamma-ray Imag- 
ing System (AGIS, Buckley et al.||2008 l, but we indeed ex- 
pect a promising future for sub- TeV detections at these high 
redshifts. 
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